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Threshold  —.2  .......  . 

Fig.  213  Variation  on  Boolpass  Operation  1. 

Image  Field  is  Replaced  (as  Input  TWO 
to  the  AND  Operation)  by  7x7  Low  Pass 
Filtered  Image  of  Field  (Fig.  119). 
Threshold  =  .2  . . . 

Fig.  214  Variation  on  Boolpass  Operation  1. 

11x11  Low  Pass  Filtered  Image  of  Field 
is  Replaced  (as  Input  ONE  to  the  AND 
Operation  )  by  a  7x7  LPF  Image  of  Field 
(Fig.  119).  Threshold  =.5  . 

Fig.  215  Variation  on  Boolpass  Operation  1. 

11x11  Low  Pass  Filtered  Image  of  Field2 
is  Replaced  (as  Input  ONE  to  the  AND 
Operation)  by  a  7x7  LPF  Image  of  Field2 
(Fig.  126).  Image  Field2  is  Replaced 
(as  Input  TWO  to  the  AND  Operation)  by 
5x5  LPF  Image  of  Field2.  Threshold  =.3  .  .  . 

Fig.  216  Variation  on  Boolpass  Operation  1. 

Image  Dispersal  is  Replaced  (as  Input 

TWO  to  the  AND  Operation)  by  7x7  Low 
Pass  Filtered  Image  of  Dispersal 

(Fig.  148).  Threshold  =.2  . 

Fig.  217  Original  Digitized  Image  Truck  (Fig. 

7).  Image  Truck  is  Replaced  (as  Input 
TWO  to  the  AND  Operation)  by  7x7  Low 
Pass  Filtered  Image  of  Truck. 

Threshold  =  .4 . . 

Fig.  218  Original  Digitized  Image  Tomcat  (Fig. 

48).  Image  Tomcat  is  Replaced  (as 

Input  TWO  to  the  AND  Operation)  by  7x7 
Low  Pass  Filtered  Image  of  Tomcat. 
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Threshold  =0  . 

Same  as  Fig.  218  Except  Threshold 
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220  Same  as  Fig.  219  Except  Threshold  =  .4  .  .  .  298 

221  Multiple  ’AND'  Operation  on  Dispersal, 

Fig.  73.  Threshold  =  .3 . 299 

222  Original  Digitized  Image  -  Field 

(Fig.  53).  Boolpass  Operation  1  Output 
Image,  Fig.  213  OR’ed  Wih  Boolpass 
Operation  2  Output  Image,  Fig.  124  .....  300 

223  Original  Digitized  Image  -  Dispersal 

(Fig.  73).  Boolpass  Operation  1  Output 
Image,  Fig.  76  OR’ed  With  Boolpas 
Deration  2  Output  Image,  Fig.  153  .....  301 

Tank  (Fig.  68)  Processed  by  Boolpass 
Operation  3 
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Abstract 


/ 

i  ' 

's  Several  new  and  seemingly  successful  scene  analysis 
techniques  for  application  to  'real’  image  processing  are 
presented  and  discussed.  These  techniques  consist  of 
particular  combinations  of  spatial  low  pass  filtering, 
global  thresholding  and  Boolean  operators,  specifically  the 
'AND’,  'OR’  and  'NOT’  operators.  These  combinatorial 
operators,  hereafter  referred  to  as  the  "Boolpass" 
operators,  perform  the  task  of  picture  ener gy/i nf ormat i on 
reduction,  while  retaining  the  fundamental  picture 

primitives  such  as  edges  whjch  characterize  the  images. 

< 

Over  150  figures  are  included  which  illustrate  the 
results  obtained  from  application  of  the  Boolpass  technique 
to  eight  different  natural  scenes.  These  results  indicate 
that  the  Boolpass  operators  do  display  great  potential  as 
important  components  of  a  larger  more  comprehensive  pattern 
recognition  machine.  Such  a  machine  would  encompass  further 
processing  <  for  target  classification/recognition  )  of  the 
resulting  Boolpass  operator  information.  /\ 

7 

/ 
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COMBINED  SPATIAL  FILTERING  AND 


BOOLEAN  OPERATORS  APPLIED  TO  THE 
PROCESSING  OF  REAL  IMAGES 

I.  INIRODUCIION 

Over  the  past  two  decades  there  has  been  an  ever 
increasing  need  to  design  and  build  machines  capable  of 
performing  tasks  which  until  now  have  been  performed  only  by 
humans.  With  the  large  and  ever  increasing  demands  being 
placed  on  mechanization  and  information  processing, 
automation  of  such  processes  has  become  an  intense  and 
challenging  area  of  current  research. 

Attempts  at  modelling  human  functions  and  intelligence 
have  lead  to  many  new  and  diverse  fields  of  study.  One  such 
field  is  Automatic  Pattern  Recognition  and  Image/Pi cture 
Processing.  The  importance  and  emphasis  that  has  been 
placed  on  this  field  of  study  has  been  generated  because, 
potentially,  it  has  the  capability  of  revolutionizing  the 
nature  of  information  processing  in  large  areas  of 
scientific,  industrial,  social,  and  military  activities 
[Ref.  43.  At  present  pattern  recognition  and  image 
processing  have  been  utilized  in  such  areas  as: 
Commerci al /Government  -  for  automatic  detection  of  flaws, 
impurities,  etc.|  natural  resource  identification!  object 
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recognition  and  parts  handling;  character  recognition;  and 
automatic  processing  of  documents; 

Medical  Applications  -  identification,  detect i on/di agnosi s 
in  X-Ray  photographs,  el ectorcardi ogram  and 
el ectroencephal ogr am  waveforms;  and 

Military  applications  -  target  detection;  i nterpretati on  of 
aerial  reconnaissance  imagery;  guidance  applications  [Ref. 
61 

Further  examples  could  be  cited.  However,  those  given  should 
clearly  indicate  the  many  diverse  and  important  application 
areas  of  pattern  recognition  research.  Advances  in  computer 
technology  have  only  recently  provided  the  more  powerful 
hardware  and  software  tools  that  will  be  required  to  find 
possible  solutions  to  the  automatic  pattern  recognition 
probl em. 

One  of  the  major  and  most  difficult  fields  of 
contemporary  pattern  recognition  and  scene  analysis  research 
is  in  the  area  of  real  image  processing.  ’Real  Image’  is 
used  here  to  mean  real-world  scenes  as  opposed  to  modelled 
scenes  or  laboratory  manufactured  scenes  which  are  quasi 
’representative’  of  the  real  world.  Such  modelled  or 
representative  images  have  unfortunately  proven  in  the  past 
to  often  oversimplify  the  pattern  recognition  problem, 
leading  to  algorithms  with  limited  or  no  capabilities  when 
applied  to  the  'real  world'  [Ref.  13. 


Problem  and  Scoge 

Scene  analysis  and  pattern  recognition  techniques  -for 
real  world  object  recognition  have  received  a  great  deal  of 
research  and  development  effort  with  varying  degrees  of 
success  CRef.  1,  2,  4,  8,  9,  111.  A  primary  objective  of 
scene  analysis  has  been  the  development  of  methods  to  reduce 
the  very  large  quantity  of  picture  data  to  suitably  relevent 
feature  information.  These  relevent  features  may  consist, 
for  example,  of  line  segments,  edges,  regions  or  other 
features,  any  or  all  of  which  can  be  used  to  represent  the 
scene  structure  CRef  ID.  The  purpose  of  this  thesis  was  to 
use  a  hierarchial  approach  to  the  investigation  and 
development  of  digital  operators  which  would  facilitate  the 
extraction  of  relevent  two-dimensional  scene  structures  from 
digitized  real  images.  The  operators  had  to  demonstrate 
picture  independence  in  their  application  and  no  a  priori 
knowledge  of  the  picture  could  be  assumed. 

"In  image  pattern  recognition  one  of  the  most  difficult 
steps  is  to  extract  objects  from  an  irrelevant  background. 
The  degree  of  difficulty  involved  varies  greatly  with  the 
quality  of  the  picture  and  the  nature  of  the  object."  CRef 
3D 

In  this  thesis,  several  new  and  seemingly  successful 
scene  analysis  techniques,  for  application  to  'real’  image 
processing,  are  presented  and  discussed.  These  techniques 
consist  of  particular  combinations  of  spatial  low  pass 


■filtering. 


gl obal 


thresholding  and  Boolean  operators. 


specifically  the  'AND',  'OR'  and  'NOT'  operators.  These 
combinatorial  operators  were  specifically  chosen  so  that 
their  simplicity  and  speed  of  operation  would  allow  their 
adaptation  to  a  real  time  processing  environment  via 
hardware  implementation.  These  combinatorial  operators, 
hereafter  referred  to  as  the  "  Boolpass  ”  operators 
(explained  in  Section  III),  are  not  to  be  interpreted  as 
independent  pattern  recognition  operators. 

The  task  of  the  Boolpass  operators  is  one  of  picture 
information  reduction  in  real  images  while  retaining  the 
fundamental  picture  primatives  such  as  edges  which 
characterize  the  images.  The  Boolpass  operators  da  display 
great  potential  as  important  components  of  a  larger  more 
comprehensive  pattern  recognition  machine.  Such  a  machine 
would  encompass  further  processing  of  the  resulting  Boolpass 
operator  information.  Further  processing  could  include 
template  matching,  moment  analysis,  size  and  range 
d i scr i mi nat i on  operators,  or  other  relevent  operators  in 
order  to  achieve  object (s)  recognition. 

General.  Approach 

In  order  for  a  computer  to  be  able  to  work  with  a 
picture  (or  image),  the  picture  must  first  be  transformed 
into  a  'digital'  image  and  placed  on  a  suitable  storage 
medium.  This  task  was  accomplished  in  the  Air  Force 
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Institute  of  Technology  <AFIT)  Signal  Processing  Laboratory 
by  use  o-f  a  video  digitizer.  The  digitizer  system  consisted 
of  a  vidicon  television  camera,  an  analog  to  digital 
converter  to  convert  the  signals  obtained  -from  the  vidicon 
camera,  and  a  magnetic  disk  storage  unit  to  store  the 
resulting  digital  picture  array  for  subsequent  processing. 
The  array  representing  the  digitized  picture  consisted  of 
256  by  256  ’pixel’  elements.  Each  pixel  element  was  a  number 
from  0  to  15  which  represented  the  grey  level  or  brightness 
of  a  very  small  region  in  the  scanned  image.  In  this  manner 
each  pixel  was  represented  in  the  computer  by  a  four  bit 
integer  value. 

The  video  digitizer  was  used  to  digitize  a  selection  of 
real  scene  images  which  were  then  stored  on  disk  for  further 
processing.  Significant  additional  noise  was  introduced  into 
the  digitized  video  images  from  both  the  A/D  converter  and 
the  camera  itself.  This  additional  noise  was  characterized 
by  inconsistent  or  spurious  pixel  value  information  randomly 
distributed  throughout  the  digitized  image.  Ambient  light 
conditions  also  affected  the  quality  of  the  digitized  image. 
In  order  to  improve  the  signal  to  noise  ratio  of  the 
digitized  picture,  a  series  of  seven  di gi ti rati ons  of  each 
picture  were  made  for  each  photograph  being  analyzed.  These 
seven  pictures  were  then  averaged  ,  pixel  by  pixel  over  the 
seven  pictures,  and  the  resultant  image  stored  as  the 


initial  image  for  processing. 


The  overall  development  of  the  'Bool  pass  operators’ 
consisted  of:  (1)  the  selection  of  a  simple  low  pass 
spatial  filter  which  could  be  applied  ’locally’,  (as  opposed 
to  a  global  filter);  <2>  the  selection  of  filter  or  ’mask’ 
sizes  to  be  applied  to  the  picture;  (3)  a  method  of 
thresholding;  (4)  the  retention  of  the  original  pixel 
values  throughout  the  processing  (  as  opposed  to  a  binary 
’black  and  white’  representati on  after  processing);  (5)  and 
an  iterative  approach  to  examining  the  potential  of  various 
combinations  of  low  pass  filtering  using  Boolean  AND 
operators  and  of  particular  importance  -  both  POSITIVE  and 
NEGATIVE  images. 

A  series  of  interactive  computer  programs  were 
developed  for  processing  the  digitized  images.  These 
programs  allowed  the  user  to  very  quickly  and  easily 
manipulate  the  manner  and  parameters  by  which  a  picture  was 
being  processed.  In  this  way  the  hierarchical  approach  to 
scene  analysis  was  accelerated  and  the  use  of  interactive 
programs  proved  to  be  a  significant  contributing  factor  in 
the  development  of  the  work  reported  in  this  thesis.  A 
detailed  description  of  the  function  of  the  computer 
programs  used,  along  with  a  Pascal  listing  for  each,  can  be 
found  in  Appendix  D. 

S§9yence  of  Present  at  i_on 

Section  II  of  this  thesis  provides  a  brief  review  and 
discussion  of  some  of  the  major  sources  of  noise  and 
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distortion  which  occur  in  a  digitized  image  as  a  result  of 
the  digitization  process.  The  purpose  o-f  this  discussion  is 
to  provide  the  reader  with  a  better  understandi ng  of  some  of 
the  complexities  and  difficulties  encountered  when  working 
with  digitized  images,  and  in  particular  'natural'  or  'real' 
digitized  images. 

A  detailed  explanation  of  the  type  and  functioning  of 
the  operators  used  in  the  thesis  is  given  in  Section  III. 
The  operators  described  in  this  section  include  low  pass 
filtering,  the  Boolpass  operators,  and  thresholding.  The 
results  obtained  from  processing  the  digitized  pictures  by 
applying  specific  combinations  of  these  operators  appear  in 
Section  IV,  Appendix  A,  B  and  C.  Both  aerial  and  ground 
scenes  were  considered  and  used  for  the  operator  testing, 
and  examples  of  each  appear  in  these  sections. 

Conclusions  and  recommendations  for  future  research  are 
discussed  in  Section  V.  Appendix  D  contains  the  listings 
and  documentation  for  the  interactive  computer  programs 
used . 
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II.  PROBLEMS  ASSOCIATED  WITH  REAL  DIGITAL  IMAGE  PROCESSING 

Noise  And  Distortion  In  Digitized  Pictures 

Noise  in  digitized  images  produces  variations  in  both 
the  intensity  and  position  of  the  pixel  points.  This  noise, 
which  is  normally  considered  undesirable,  is  introduced  into 
the  digitized  picture  from  several  sources. 

The  first  source  to  be  considered  is  the  picture, 
itself,  which  is  being  digitized.  Noise  can  already  exist  in 
a  picture  image  or  photograph  due  to  poor  picture  quality, 
processsing  or  reproduction.  Improper  lighting  can  produce 
distortion  and  improper  film  speed  can  produce  a  grainy’  or 
blurred  effect  on  a  photograph.  These  and  other  effects 
introduce  ’noise’  into  the  picture  and  detract  from  the 
ideal:  a  clear,  crisp  image  desired  for  processing. 

The  digitizer  is  most  frequently  the  major  source  of 
noise  in  the  image  processing  system.  The  digitizer 
consists  of  five  main  elements:  (1)  the  sampling 
aperature,  which  accesses  the  individual  pixel  areas;  (2) 
the  scanning  mechanism,  which  sequentially  moves  the 
sampling  aperature  over  the  image;  (3)  the  sensor,  which 
measures  the  brightness  of  each  pixel  (a  linear  response  is 
highly  desirable  here);  (4)  an  analog  to  digital  converter 
to  quantize  the  output  of  the  sensor;  and  (5)  an  output 
medium  for  storage  and/or  display.  (If  required,  a  detailed 
description  and  discussion  of  each  can  be  found  in  Ref.  2.) 
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Each  of  the  above  elements  in  a  digitizer  contributes 
to  the  amount  of  noise  that  will  be  added  to  the  digitized 
image.  Total  elimination  of  this  noise  is  not  physically 
possible.  Therefore  an  image  initially  digitized  by  such  a 
system  is  inevitably  degraded  before  the  actual  processing 
even  begins. 

"If  a  uniformly  gray  image  is  presented  to  a  digitizer, 
the  noise  inherent  in  the  system  will  cause  variations  in 
the  output  even  though  the  input  is  constant.  Noise 
introduced  by  the  digitizer  is  a  source  of  image  degradation 
and  should  be  small  relative  to  the  contrast  of  the  image." 
CRef.  21 

In  natural  or  real  scenes  the  contrast  of  the  image  and 
the  amount  of  energy  in  the  image  can  vary  widely.  In  such 
pictures,  image  processing  becomes  an  extremely  difficult 
task  as  the  additional  presence  of  noise  can  effectively 
hide  or  distort  relevent  boundary,  edge,  and  texture 
information  which  was  present  in  the  original  scene.  It  is 
believed  that  these  fundamental  primatives  -  boundary,  edge, 
texture  and  color  information  are  possible  keys  to  the 
solution  of  the  pattern  recognition  problem  LRef.  93. 
Finding  operators  which  can  be  applied  universally  to 
natural  images,  where  such  key  features  as  micro  and  macro 
edges  are  possibly  lost  or  distorted,  creates  one  of  the 
most  complex  and  as  yet  unsolved  problem  areas  in  image 
processing.  (  A  micro  edge,  for  example,  could  be  defined 
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as  two  adjacent  points  having  different  grey  levels,  while  a 
macro  edge  might  exist  between  two  large  coarsely  textured 
regions  having  different  average  pixel  greylevels  CRef  103.) 

Examples  of  digitised  pictures  using  the  AFIT  Vidicon 
system  are  shown  in  figures  1  thru  3.  In  these  examples, 
the  existence  of  edge  noise,  distortion  and  discontinuities 
are  clearly  visible.  However,  the  fact  that  the  human 
visual  system  can  still  readily  discern  distinct  objects  in 
these  pictures,  despite  this  noise,  leads  one  to  believe 
that  object  extraction  by  a  machine  with  this  same  picture 
information  should  still  be  possible. 


Fig.  1. 


Digitized  Image  of  Black  Blocks  on  a  White 
Backgroud  (Note  ’Noise’  and  Distortion  Introduced 
By  Digitization  Process) 


III.  OPERATOR  DESCRIPTIONS 


Low  Pass  Fi_l_t  er 

A  noise  smoothing  procedure  often  applied  to  a 
digitized  picture,  before  feature  or  edge  extraction,  is  low 
pass  filtering  CRef.  13.  This  type  of  filtering  is 
particularly  useful  when  applied  to  real  images  where  the 
edges  frequently  suffer  from  noise  or  some  source  of 
discontinuity.  Examples  of  this  were  previously  illustrated 
in  Section  II  (figures  1  thru  3). 

Low  Pass  Filtering  directly  affects  the  gradient  of  a 
picture  by  smoothing  out  fast  gradients  and  reinforcing  slow 
gradients.  The  gradient  magnitude  operator  basically 
represents  the  steepness  of  slope  at  every  point.  In  areas 
of  steep  slope,  such  as  at  edges,  the  gradient  takes  on 
large  values,  while  low  gradient  values  result  in  areas  of 
nearly  uniform  grey  level  distribution.  Application  of  the 
low  pass  filter  results  in  a  blurred  image  and  suppression 
of  isolated  point  noise  sources. 

Two  low  pass  filters  were  used  in  this  thesis.  Each 
filter  was  computed  by  a  simple  averaging  of  the  image 
points  for  a  given  mask  size.  The  filter  ’mask’  refers  to 
the  matrix  of  pixels  being  processed  by  the  filter. 

The  first  of  the  low  pass  filters  used  was  a  two 
di mensi onal (2-D)  filter  which  operated  in  both  the  vertical 
and  horizontal  axis  of  the  picture  plane. 

1 A 


The  second  low 


pass  -filter  was  one  -  di  mensi  onal  <  1  -D>  and  operated  only  in 

the  horizontal  axis  of  the  picture.  Five  different  mask 

sizes  were  used  and  tested  for  each  filter,  namely  -  3,  5, 

7,  9  and  11  pixels.  These  numbers,  for  the  two-dimensional 

filter,  represent  the  number  of  pixels  along  one  edge  of  a 

square  matrix.  As  an  example,  consider  image  points 

arranged  as  follows  for  a  mask  of  size  3  (or  3x3  mask): 

ABC 
D  X  E 
F  G  H 

The  result  of  the  two-dimensional  Low  Pass  Filter  (LPF) 
averaging  operation  at  point  X  would  be: 

LPF  at  X  =  (A+B+C+D+X+E+F+G+H)/9 
Using  this  same  example  for  the  one-dimensional  low 
pass  filter  for  a  mask  size  of  3  (or  1x3  mask): 

LPF  at  X  =  (D  +  X  +  E)/3 

It  can  be  seen  that  the  mask  size  specified  for  a  one- 
dimensional  low  pass  filter  represents  an  array  symmetric 
about  the  point  being  processed  (eg.  X) .  The  LPF  at  pixel  X 
(the  matrix  center)  for  the  remaining  matrix  sizes  -  5,  7, 

9  and  11,  were  computed  in  a  similar  manner. 

In  this  thesis,  the  Low  Pass  Filter  was  applied  in  the 
image  domain.  That  is,  the  LPF  operator  was  applied  to  each 
pixel  of  the  original  digitized  image  (i.e.  the  image 
domain)  as  opposed  to  the  transform  domain.  The  transform 
domain  results  when  one  or  more  operations  are  performed  on 
the  original  image,  producing  an  already  processed  or 
’transformed’  image. 
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Many  variations  and  more  complex  Low  Pass  Filters  are 
possible.  However,  the  LPF  explained  -for  use  in  this 
thesis  was  chosen  -for  its  ease  of  implementation  and  speed 
of  operation  on  almost  any  computer  system.  Examples  of  the 
filter  as  applied  to  real  images  can  be  found  in  Section  IV 
RESULTS. 

I(TC§shgl_di_ng 

Thresholding  was  used  in  conjunction  with  the  Boolpass 
operators  to  specify  the  acceptance/rejection  limits  for  any 
one  particular  operation  during  the  processing  of  the  image. 
The  thresholding  procedure  used  here  with  the  Boolpass 
operators  permitted  retention  only  of  those  gray  level 
values  within  a  certain  window  of  values,  determined  for 
each  individual  pixel  in  the  image  being  processed.  All 
other  values  were  set  to  grey  level  15  which  equals  a  blank 
(white) . 

Threshold  values  were  based  on  a  calculation  of  the 
global  mean  and  standard  deviation  of  the  pixel  values  for 
the  entire  digitized  picture. 

That  is: 

Mean  pixel  value  =  sum  of  pixel  grey  level  values  for 

entire  picture  divided  by  the  total 
number  of  pixels 

Standard  Deviation  =  Square  root  of  the  average  of  the 

squared  deviations  from  the  mean 

A  percentage  of  the  standard  deviation  was  then  used  to 


1  b 


determine  the  minimum  and  maximum  threshold  limits. 
(Mathematical  formulae  for  these  common  mean  and  standard 
deviation  computations  are  available  in  the  CRC  Standard 
Mathematical  Tables<  24th  Edition),  page  473.) 

The  choice  of  using  a  gl obal  approach  to  deter mi ne  the 
limits  for  thresholding,  as  opposed  to  a  local  or  locally 
adaptive  threshold  technique  CRef.  91  was  initially  an 
arbitrary  decision  based  on  computational  simplicity  and 
speed  of  operation.  However,  once  the  threshold  limits  were 
calculated  they  were  applied  to  each  individual  pixel  so  as 
to  determine  a  threshold  window,  or  window  of  acceptance, 
for  that  particular  pixel  (grey  level).  Subsequent 
processing  using  this  combination  of  global  threshold  limits 
together  with  a  locally  determined  threshold  window  for 
each  pixel  gray  level,  yielded  positive  results  while 
requiring  only  simple,  fast  calculations.  The  following 
example  will  help  illustrate  how,  within  the  interactive 
programs,  the  threshold  limits  and  windows  were  calculated. 

When  an  image  was  to  be  processed,  the  user  was 
prompted  at  the  computer  terminal  to  enter  a  decimal  number 
between  0  and  1.  The  standard  deviation,  having  previously 
been  computed  for  the  picture,  was  then  multiplied  by  this 
decimal  number.  The  result  of  this  multiplication  was  then 
rounded  off  to  an  integer  value.  Plus  or  minus  this  integer 
value  then  gave  the  threshold  limits  (or  global  threshold 
limits).  This  calculation  was  performed  only  once,  at  the 
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start  of  processing  of  each  picture. 


The  threshold  limits  were  then  used  at  each  pixel 
processed  in  the  image  to  calculate  a  threshold  range  or 
threshold  window  for  that  particular  pixel.  Very  simply, 
the  threshold  limits  (range  or  window)  at  each  pixel  was 
determined  by  addition  and  subtraction  of  the  (global) 
threshold  limit  to  the  gray  level  value  of  the  pixel. 

For  example: 

given  a  digitized  picture  256  >;  256  pixels  (16 

greyl evel s) 

mean  of  pixel  values  =  8 
standard  deviation  =  3 
decimal  number  entered  by  user  =  .4 

THRESHOLD  LIMITS  =  -  round (.4  x  3) ,  +  round (.4  X  3) 

THRESHOLD  WINDOW  at  pixel  P  =  P  -  round (.4  x  3)  , 

P  +  round(.4  x  3) 

Sol vi ng : 

THRESHOLD  WINDOW  (min,  max)  at  pixel  P  =P-1,  P+1 
if  P  =  10,  then 

THRESHOLD  WINDOW  (min,  max)  =9,  11 

Therefore,  if  a  Boolpass  operation  is  performed  using 
the  above  threshold  window,  only  pixel  values  of  9,  10,  and 
11,  at  pixel  P  =  10,  will  be  retained  (or  otherwise  operated 
upon)  while  all  other  values  wi 1 1  be  set  to  gray  level  15. 

The  coefficient  chosen  to  establish  the  threshold 
limits  has  a  significant  impact  on  the  results  which  are 
received  from  applying  the  Boolpass  operators.  Examples  of 
various  threshold  limits  used  with  the  Boolpass  operators 
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are  demonstrated  in  Section  IV 


r esul ts . 


Booigass  Operators  and  Thei^r  Development 

"Boolpass"  was  chosen  as  an  acronym  for  the 
combinatorial  operators  or  picture  processing  techniques 
that  evolved  from  the  research  and  testing  in  this  thesis. 
"Boolpass"  was  based  on  a  contraction  of  the  two  principle 
operations  used  for  the  image  processing  -  namely,  BOOLean 
operations  and  Low  PASS  filtering  operations. 

The  Boolpass  operators,  as  used  in  this  thesis,  were 
designed  to  operate  simultaneously  on  two  ’input’  images  and 
produce  a  third  ’output’  or  processed  image.  In  this  manner 
the  contents  (pixel  values)  of  the  input  image  files  were 
not  changed  after  a  Boolpass  operation  therefore  permitting 
additional,  separate  processing  of  the  original  image(s),  if 
requi red . 

The  Boolean  operations  used  in  this  thesis  consisted  of 
the  "AND",  "OR"  and  "NOT"  operations.  Though  frequently 
considered  in  the  context  of  binary  operations,  these 
Boolean  operators  were  not  limited  to  just  binary 
applications  here.  Instead,  Boolean  operations  were  applied 
to  any  or  all  16  pixel  gray  levels.  The  range  of 
application  of  these  Boolean  operations  was  dependent  on  (1) 
the  actual  picture  greylevel  content,  and  (2)  the  gray  level 
values  specified  by  the  threshold  window. 

A  basic  illustration  of  the  functioning  of  these 
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operations,  using  only  binary  values  ’0’  (TRUE)  and  ’1’ 
(FALSE),  is  given  by  the  -following  truth  tables.  As  was 
already  noted,  in  this  thesis  the  -following  operations  were 
extended  to  as  many  as  16  values  -  0  thru  16.  ("A  truth 
table  gives  the  truth  values  of  a  logical  expression  for 
each  combination  of  the  truth  values  of  its  variables." 
[Ref.  71 
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These  Boolean  operators  were  used  to  perform  two 
important  processing  operations:  (1)  pixel  by  pixel 
comparison  of  two  digitized  images,  and  (2)  creation  of 
inverse  greylevel  images,  ie.  negative  pictures  from 
positive  pictures  and  positive  pictures  from  negative 
pictures. 

The  first  operation,  pixel  comparison,  could  be 
performed  with  or  without  thresholding.  Without 
thresholding,  the  pixel  comparison  was  accomplished  by  only 
the  ’AND’  operator.  Using  the  AND  operation,  correspondi ng 
pixel  positions  in  the  two  input  images  were  sequentially 
compared  or  ’ANDED’.  If  the  two  values  being  compared  were 
the  same,  that  value  was  stored  in  a  third  or  output  picture 
file.  If  the  values  being  compared  were  not  the  same,  a 
value  of  15  (blank)  was  stored  in  the  output  picture  file. 
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(Note  that  the  pixel  values  that  were  stored  in  the  output 
•file  did  not  (necessarily)  result  in  a  binary  picture  -  with 
values  0  (black)  and  15  (white)  only.)  This  ANDING 
operation  was  performed  on  all  64k  picture  pixels  and  the 
results  stored  as  a  new  picture  file,  leaving  the  original 
input  images  unchanged. 

When  pixel  co.nparison  was  performed  with  thresholding, 
both  the  ’AND'  and  ’OR’  operators  were  required.  The  OR 
operation  was  needed  because  more  than  one  value,  as 
specified  by-  the  threshold  window,  was  acceptable  during  the 
AND  operation.  In  performing  the  pixel  comparison  with 
thresholding,  a  single  pixel  value  in  one  input  image  file 
was  compared  to  a  range  of  pixel  values  (as  specified  by  the 
threshold  window)  for  the  corresponding  pixel  in  the  second 
input  image  file.  The  OR  operator  was  then  used  for 
comparison  of  this  range  of  values  in  conjunction  with  the 
AND  operation.  Further  discussions  using  this  combination 
of  AND  and  OR  operations  will  be  simply  referred  to  as  "AND" 
or  "ANDING".  The  existence  of  the  OR  operation  will  be 
implicit  to  any  "ANDING"  operation  where  Threshold  imposes  a 
range  rather  than  a  single  value  for  the  comparison.  The 
following  example  should  help  illustrate  these  pixel 
comparison  concepts. 

Consider  two  input  file  images  -  Ain,  Bin,  and  an 
output  file  image,  Pout.  Corresponding  pixel  locations  in 
each  are  represented  by  A1 ,  B1  and  PI.  For  this  example  let 


A1  =  6  and  Bl  =  7.  Threshold  limits  were  calculated  at  .5 
(interactive  user  input)  X  standard  deviation  to  be  equal  to 
the  pixel  value  plus  or  minus  2.  The  ANDING  operation  is 
now  performed  between  A1  and  Bl: 

First,  without  thresholding  :  given  A1  =  6  and  Bl  =  7, 
does  A1  =  6  AND  Bl  =  6?  Obviously  not,  therefore  a  value  of 
15  is  stored  in  PI. 

Next,  with  thresholding  :  given  A1  =  6,  Bl  =  7  and  a 
threshold  limit  of  +  or  -  2,  does 

A1  =6  AND  Bl  =  4 
OR 

A1  =6  AND  Bl  =  5 
OR 

A1  =6  AND  Bl  =  6 
OR 

A1  =6  AND  Bl  =  7 
OR 

A1  =6  AND  Bl  =  8 

Since  A1  =  6  AND  Bl  =  6,  for  the  threshold  window  given 
at  Bl,  the  answer  is  yes,  and  so  a  pixel  value  must  be 
selected  and  stored  in  PI.  The  value  of  the  pixel  chosen  to 
be  stored  in  Pout,  for  example  -  either  A1  or  one  of  the 
five  threshold  window  ’B’  values,  can  have  a  significant 
impact  on  subsequent  processing  of  Pout.  Details  of  this 
type  of  processing  will  be  discussed  later  in  this  section. 

The  initial  research  of  this  thesis,  into  the  problem 
of  processing  real  digitized  images,  concentrate,-  on  a 
search  for  a  universally  applicable  method  which  could  be 
used  to  help  reduce,  or  at  least  more  evenly  distribute,  the 
large  amount  of  noise  that  was  always  present  in  the  picture 
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after  digitization.  Low  pass  filtering,  as  previously 
described,  was  found  to  be  very  useful  in  helping  to  solve 
this  noise  distribution  problem. 

The  filter  or  mask  size  used  in  the  low  pass  filter 
<LPF)  operation,  during  picture  processing,  resulted  in  very 
different  output  picture  images.  Examples  of  these 
differences  using  3x3,  7x7,  and  11x11  mask  sizes  are 
illustrated  below  in  figures  4  thru  6). 
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During  the  research,  examination  of  these  and  many 
other  examples  showed  a  definite  tendency  of  object  or 
feature  edge  structures  to  remain  while  more  inconsistently 
textured,  noisy  areas  were  smoothed.  The  discovery  that 
certain  basic  structures  remained  relatively  unchanged  after 
processing,  while  others  did  not,  provoked  interest  over 
what  features  did  remain  unchanged  after  processing  with 
different  size  masks.  It  was  thought  that  a  direct 
comparison  of  two  low  pass  filtered  images  might  provide  an 
answer  to  the  question. 

The  images  used  for  comparison  were  identical  except 
that  each  had  been  previously  processed  by  a  different  size, 
but  same  dimension,  low  pass  filter.  Comparison  of  the  two 
images  was  accomplished  through  use  of  the  Boolean  AND  (and 
OR)  operation.  The  ANDING  operation,  depending  on  the 
threshold  setting,  was  able  to  substantially  reduce  the 
amount  of  energy  in  the  picture.  However,  feature  extraction 
at  this  level  was  still  considered  difficult  and  further 
testing  using  combinations  of  different  size  filtered 
images,  thresholding,  and  ANDING  gave  interesting  but  mixed 
results.  For  examples  see  Section  IV  -  Results. 

Large  areas  af  unwanted  picture  information  still 
remained  in  those  areas  where  the  picture  gradient  was  very 
small.  Attempts  to  reduce  these  areas  using  picture 
thresholding  and  conventional  edge  operators  (such  as  Kirsch 
and  5-level  simple  masks  CRef.  93),  met  with  no  success. 
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Further  visual  study  of  the  digitized  pictures  showed  that 
a\_  least  some  of  this  unwanted  picture  energy  might  be 
removed  if  the  undesired  area,  of  a  particular  greylevel, 
was  combined  with  its  complement.  The  solution  to  this 
problem  was  the  creation  and  use  of  negative  images  during 
the  combinational  image  comparisons. 

For  any  given  image,  the  complement  or  negative  of  that 
image  ANDED  to  itself  would,  of  course,  yield  an  output 
image  consisting  only  of  blanks.  The  problem  which  then  had 
to  be  solved  was  what  negative  image  could  be  used  to  help 
cancel  or  remove  the  ’undesired'  areas  in  the  image  being 
processed.  The  solution  to  this  problem  was  to  use  a  low 
pass  filtered  image  of  a  different  mask  sire  than  that  of 
the  image  being  processed.  These  filtered  images  contained 
similar,  but  not  exactly  the  same  pixel  information  to  that 
of  the  image  being  processed.  The  negative  of  this  filtered 
image  was  then  used  as  one  of  the  input  images  to  the 
Boolean  AND  operation.  The  original  image  that  was  being 
processed  was  used  as  the  second  input  image  to  the  AND 
operation.  Therefore,  when  a  positive  image  and  a  low  pass 
filtered  image  were  simultaneously  processed  by  the  AND 
operation,  a  large  reduction  in  pixel  energy  became 
possible.  The  amount  of  energy  reduction  achieved  was 
dependent  on  the  size  and  type  of  filter  used,  as  well  as 
the  threshold  window  chosen  for  the  Boolean  operation. 

Application  and  experimentation  with  this  technique  of 
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ANDing  selective  positive  and  negative  images  proved  to 
yield  some  very  interesting  results.  Test  results  showed  a 
significant  improvement  in  total  picture  energy  reduction 
while  relevent  feature  information  was  retained.  Detailed 
results  are  illustrated  in  Section  IV. 

It  has  already  been  shown,  though  a  previous  example, 
that  during  the  processing  of  two  images  by  the  AMD 
operator,  pixel  values  are  stored  in  a  separate  picture  file 
to  form  the  new  processed  picture.  During  the  combinational 
testing  of  the  Boolpass  operators,  the  values  which  were 
chosen  to  be  stored  and  represent  the  new  processed  picture 
were  significant  to  any  subsequent  processing  done  to  that 
picture.  The  reason  for  this  can  be  attributed  to  the 
threshold  window  specified  for  the  subsequent  processing 
operations.  The  following  example  should  best  illustrate 
why. 

For  this  example,  consider  three  input  image  files  — 
Ain,  Bin,  Cin,  and  three  corresponding  pixel  locations  in 
each  -  Al,  Bl,  Cl.  If  A1  was  being  ANDED  with  B1  and  the 
values  were  exactly  equal,  it  would  not  matter  from  which 
image  the  pixel  value  was  taken  to  be  stored  into  the 
processed  picture  matrix.  Pout.  However,  if  the  pixel 
values  were  not  exactly  equal,  acceptance  or  rejection  is 
based  upon  the  limits  specified  by  Threshold.  Again,  let  Al 
=6,  Bl  *  7,  and  Cl  =  9.  Threshold  limits  were  calculated 
at  .5  X  standard  deviation  to  be  equal  to  the  pixel  value 
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plus  or  minus  2. 

The  ANDING  operation  is  now  performed  between  A1  and 
Bl,  ie  given  A1  =  6  (and  a  threshold  limit  of  +  or  -  2), 
does  Bl  =  4  OR  5  OR  6  OR  7  OR  B  ?  Since  Bl  =  7,  the  answer 
is  yes,  and  so  a  pixel  value  must  be  selected  and  stored  in 
’Pout’.  I-f  the  value  o-f  A1  is  chosen,  PI  =  6;  if  Bl  is 
chosen,  PI  =  7.  Many  other  methods  of  pixel  value  selection 
could  be  used  based  on  any  number  of  different  criterea.  In 
this  thesis,  if  the  result  of  the  AND  operation  required 
storage  of  some  pixel  value  from  the  Bin  file,  the  Boolpass 
operators  always  stored  the  pixel  value  of  the  original 
input  ’ B ’  file.  The  choice  of  the  ’B’  input  pixels  was 
based  an  hierarchical  test  results  obtained  during 
interactive  tests  using  several  ad  hoc  criterea. 

Next,  the  same  processing  operation  is  to  be  applied 
to  a  multiprocessing  of  the  pictures  ’Cin’  and  ’P out’.  From 
the  value  given  for  Cl  (ie.  9) ,  performing  the  AND  operation 
exceeds  the  Threshold  limits  if  A1  had  previously  been  the 
value  for  retention  and  storage  in  Pout.  Therefore,  the 
value  stored  in  another  new  picture  matrix,  say  P’ ,  would  be 
a  blank.  However,  if  Bl  had  been  stored  in  Pout,  the  value 
would  have  been  within  the  threshold  limits  and  a  greylevel 
other  than  a  blank  would  have  been  stored. 

From  this  example,  it  is  clear  that  selection  criterea 
for  pixel  values  stored  after  processing  does  affect 
subsequent  processing  of  that  image.  It  is  also  for  this 
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reason  that  multi-level  processing  o-f  pure  binary  images 
using  the  Boolpass  operators  is  not  very  effective.  This  is 
due  to  the  large  amount  of  feature  information  which  is  lost 
for  subsequent  processing  as  a  result  of  reducing  the 
picture  to  a  binary  image. 
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Ith  RESULTS 


Gener  al_ 

The  number  of  combinations  of  processing  techniques 
which  are  possible  using  the  concepts  and  techniques  of  the 
Bool  pass  operators  is  very  large.  The  testing  results  which 
follow  should  not  be  considered  an  exhaustive  study  of  these 
combinations;  however,  in  exercising  the  Boolpass  operators 
on  a  total  of  8  original  images,  over  400  pictures  were 
processed.  The  development  of  the  Boolpass  operators  was 
not  the  result  of  some  totally  preconceived  theoretical  or 
abstract  concept  but  rather,  of  an  hierarchical  approach 
based  on  ’ man-i n- the— 1 oop ’  real  time  interactive  imiage 
processing  using  a  digital  computer.  The  ’loop’  referred  to 
here  consisted  of  interactive  image  processing  programs  on  a 
digital  computer;  a  video  monitor  which  displayed  the 
processed  (output)  images;  and  a  ’man’  or  ’user'  who 
controlled  the  processing  environment  and  input  parameters 
to  the  computer  programs.  The  following  diagram  illustrates 
the  interactive  ’loop’  concept. 
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This  interactive  approach  to  image  processing  allowed 
the  Lis er  to  visually  analyze  the  output  images  after  each 
stage  of  the  processing.  The  visual  analysis  then  resulted 
in  a  subjective,  but  discerning  human  'feedback'  loop  within 
the  image  processing  system.  In  this  way,  changes  to  the 
input  parameters  were  made,  if  considered  necessary,  and 
used  for  subsequent  processing.of  the  images. 

A  particular  sequence  of  processing  steps,  or  Boolpass 
operation,  was  considered  noteworthy  if,  when  applied  to 


several  very  different  images,  that  same  operation(s) 
performed  well  at  extracting  feature  information  (considered 
"useful"  by  the  user)  in  all  the  test  images. 

The  results  illustrated  herein  represent  those 
combinations  of  the  Bool  pass  operators  which  appear  to  have 
the  greatest  potential  for  use  in  an  image 
processing/pattern  recognition  machine.  It  should  be  noted 
that  the  processed  pictures  in  this  section  are  the  result 
of  applying  only  the  Boolpass  operators  (as  previously 
described  in  Section  III),  unless  specified  otherwise. 
Further  noise  reduction  and  target  extraction  from  the 
processed  images  which  follow  would  be  facilitated  by  the 
use  of  additional  target  or  feature  information.  Additional 
processing  data  such  as  range  to  target  information,  target 
size  and  feature  templates  would  be  valuaLle  aides  to 
further  refinement  of  the  images  or  target  extraction. 

It  should  be  emphasized,  again,  that  the  results  of  the 
Boolpass  operators,  as  given  in  this  thesis,  represents  only 
the  first  stage  of  a  two  stage  pattern  recognition  machine. 
This,  the  first  stage,  accomplishes  overall  picture  energy 
reduction  and  feature  edge  extraction.  The  second  stage, 
target  recognition,  is  then  performed  on  the  output  image  of 
the  first  stage  of  processing.  The  advantage  to  this 
technique  is  that  the  complex  task  of  target  recognition  is 
only  applied  to  a  small  number  of  probable  target  areas  (as 
determined  by  the  Boolpass  oerators)  ,  as  opposed  to  the 
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entire  image.  This  aids  in  reducing  the  complexity  of  the 
target  recognition  task,  which  in  turn  enhances  the 
operating  speed  of  the  overall  pattern  recognition  machine. 

of  Pc?§ent at i^on  3  Results 

The  results  have  been  divided  into  three  main  topic 
areas:  1)  low  pass  filtering  (both  one  and  two 
dimensional);  2)  boolean  operators  (  AND,  OR,  NOT);  and 
3)  Boolpass  operators  (with  thresholding).  Each  picture, 
or  sequence  of  pictures  presented,  has  been  preceeded  by  an 
outline  of  those  steps  and  input  values  which  were  used  to 
process  the  input  image(s)  to  the  final,  output  picture 
shown.  Subsection  3).  on  the  Boolpass  operators,  also 
includes  three  dimensional  (3-D)  picture  energy  plots  of  the 
original  input  and  final  output  images  to  the  Boolpass 
operators.  These  3-D  plots,  which  graphically  show  the 
pixel  value  distribution  over  the  entire  picture,  help  to 
illustrate  the  effectiveness  of  the  Boolpass  operators  in 
achieving  picture  energy  reduction. 

For  ease  of  reference,  page  vi  begins  a  detailed 
listing  of  all  pictures  (figures)  displayed  in  this  section. 

II  Low  Pass  Fil_teri_ng 

Figures  7  thru  IB  show  examples  of  two  dimensional  low 
pass  filters,  using  mask  sizes  of  3x3,  5x5,  7x7,  9x9  and 
11x11.  One  dimensional  filtering  (in  the  horizontal  axis) 
is  illustrated  in  figures  19  thru  26  for  mask  sizes  3,  5,  7 
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and  1 1 . 


It  should  be  noted  that  the  amount  of  computation 
required  for  a  two  dimensional  filter  is  at  least  N  times 
that  required  for  a  one  dimensional  filter  of  size  N.  (For 
the  computer  used  in  this  thesis,  a  Data  General  Eclipse 
S/250,  this  translated  into  the  two  dimensional  filter 
requiring  almost  3  times  the  processing  time  required  for 
the  one  dimensional  filter.)  As  might  be  expected,  further, 
but  similar  processing  of  the  output  of  these  two  filters 
also  yields  two  distinct  images.  The  one  dimensional 
filter,  because  of  computation  speed,  might  seem  the  obvious 
choice  for  a  pattern  recognition  machine.  However,  the 
choice  of  which  filter  should  be  used  in  a  particular 
algorithm  or  pattern  recognition  machine  depends  on  the 
constraints  imposed  by  that  machine  or  algorithm.  For 
example,  a  target  extraction  algorithm  may  work  very 
effectively  on  the  output  of  a  two  dimensional  filter,  but 
not  as  well  on  the  one  dimensional  ouput.  In  another  case, 
speed  may  be  critical  and  the  one  dimensional  filter  must  be 
used,  even  if  more  errors  are  introduced  into  the  target 
extraction  algorithm. 

Later  in  the  results,  examples  are  given  of  Poolpass 
operators  using  both  types  of  low  pass  filters.  The  results 
show  that  subsequent  processing  using  the  Poolpass 
operators  yield  similar  ouput  images.  This  observation 
could  imply  that  the  faster,  one  dimensional  filtering  might 
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Fig.  11.  Truck  (Fig.  7)  Processed  With  LPF  9x9 


0,2 
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Truck  (Fig.  7)  Processed  With  LPF  1x3 


I 

i 

( 


Fig.  25.  Field  <Fig.  13)  Processed  With  LPF  1x7 
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2).  Boolean  Operators 

A  description  of  the  boolean  operations  AND,  OR  and 
NOT,  as  applied  to  image  processing,  was  previously  detailed 
in  Section  III.  Further  explanation,  beyond  that  already 
given,  is  not  considered  necessary  here  to  be  able  to 
understand  and  interpret  the  processed  image  results  which 
foil ow. 

Examples  of  images  processed  using  boolean  operations  , 
both  with  and  without  thresholding,  are  presented.  This 
provides  the  reader  with  an  opportunity  to  compare  and 
contrast  the  effect  of  imposing  this  additional  constraint 
during  any  particular  boolean  operation. 
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Boolean  NOT  Operation  Performed  on 


Fig.  28.  Boolean  NOT  Operation  Performed  on  Field  (Fig.  13) 
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29.  5>: 5  LPF  Truck  (Fig.  9)  ANDED  With  Truck  (Fig. 

Threshold  =  0 


61 


Fig.  30.  5x5  LPF  Truck  iFig.  9)  ANDED  With  Truck  (Fig.  7) 

Threshold  =  .2 
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Fig.  31.  5x5  LPF  Truck  (Fig.  9)  ANDED  With  Truck 

Threshold  =*  .4 
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Fig.  34.  5x5  LPF  Truck  (Fig.  9)  ANDED  With  Truck  (Fig.  7) 

Threshold  =  1.0 
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Fig. 


mm 


36.  1x5  LPF  Field  (Fig.  24)  ANDED  With  Field  (Fig.  1) 

Threshold  =  .2 
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37.  1x5  LPF  Field  (Fig.  24)  ANDED  With  Field  (Fig.  1> 

Threshold  =  .4 
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Fig.  39 . 


1x5  LPF  Field 


(Fig.  24)  ANDED  With  Field  (Fig. 
Threshold  =  .8 


1) 
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Fig.  41 


Digitized  Image  of  Black  Triangle  on  White 
Background 


►  'V 


Fig.  42.  Boolean  OR  Operation  Performed  Between  Digitized 
Blocks  (Fig.  1)  and  Digitized  Triangle  (Fig.  42) 
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32.  ?9°1bs§§  SBgC§tgrs 

The  Boolpass  operators  use  specific  combinations  of 
those  processing  operations  demonstrated  in  the  preceeding 
two  subsections.  A  detailed  explanation  of  the  development 
and  functioning  of  the  Boolpass  operators  was  already  given 
in  Section  III  and  will  not  be  reiterated  here. 

Appendix  A  contains  examples  of  processed  image  results 
using  two  different  Boolpass  operators.  These  operators, 
outlined  below  and  labelled  simply  Boolpass  operator  1  and 
Boolpass  operator  2,  are  applied  with  the  same  parameters 
and  threshold  to  a  number  of  different  images.  This 
approach  to  presenting  the  results  permits  a  comparison  of 
each  operator’s  effectiveness  when  applied  to  a  variety  of 
input  images.  The  results  presented  in  Appendix  A  include 
both  types  o f  low  pass  spatial  filtering;  one  di mensi onal 
and  two  dimensional  <  as  demonstrated  in  subsection  1)). 
Particular  note  should  he  taken  of  how  the  threshold  value 
(used  in  conjunction  with  the  mean  and  standard  deviation  of 
the  original  picture)  affects  the  final  processed  image. 

The  following  two  image  processing  (procedural) 
outlines  detail  the  exact  processing  steps  by  which  the 
Boolpass  operator  results  in  Appendix  A  were  obtained.  The 
generalized  forms  such  as  Fig.  A,  Fig.  B,  etc.  are  replaced 
by  actual  figure  references  in  the  abbreviated  outlines 
preceeding  each  set  of  results  in  Appendix  A. 
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AO-AUQ  076  AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OH  SCHOO— ETC  F/G  20/6 

COMBINED  SPATIAL  FILTERING  AND  BOOLEAN  OPERATORS  APPLIED  TO  THE— ETC(U» 
JUN  82  B  E  FELTMATE 

UNCLASSIFIED  AFIT/GCS/EE/82J-8  NL 


Qy feline  for  Bool qsss  Oger afeor  1 


1.  Original  picture  processed  :  Fig.  A 
Step  1 

Apply  2-dimensional  11x11  low  pass  filter  to  original 
picture,  Fig.  A 

2.  Resulting  output  image  :  Fig.  B 
Step  2 

Form  negative  image  of  Fig.  B  (ouput  image  from  step  1) 

3.  Resulting  ouput  image  :  Fig.  C 

Step  3 

Perform  ’AND'  operation  using  Fig.  C  as  input  ’ONE’  to 
the  AND  operator  and  Fig.  A  as  input  ’TWO’  to  the  AND 
operator 

Case  1.  Threshold  =  (1/2  the  Standard  Deviation  of 

the  pixel  values  for  the  entire  picture. 
Fig.  A)  /  10 

4.  Resulting  ouput  image  s  Fig.  D 

Case  2.  Threshold  =  (the  Standard  Deviation  of  the 

pixel  values  fot  the  entire  picture.  Fig. 
A)  /  10 

5.  Resulting  ouput  image  :  Fig.  E 

NOTE  :  Pixels  stored  in  output  image  of  AND  operation 
are  taken  from  input  ’TWO’ 

6.  Other  images  :  as  specified 


This  same  operation  is  repeated  exactly  for  a  1 -Di mensi onal 
low  pass  filter.  The  only  change  in  the  processing 

procedure  occurs  in  Step  1,  as  follows: 

Step  1 

Apply  a  1-Dimensional  lxll  low  pass  filter  to  original 
picture,  Fig.  A 


An  abbreviated 

form 

of 

thi  s 

out  1 i ne. 

listing  only  the 

nui  ibers  1  thru  6 

and 

the 

corresponding 

correct  Figures,  is 

included  before 

the 

resul ts 

for  any 

particular  image 

processed . 

Out  l_i_ne  for  Bgolgass  Operator  2 

1.  Digitized  picture  processed  :  Fig.  A 
Step  1 

Apply  2— di mensi anal  7x7  low  pass  filter  to  original 
picture.  Fig.  A 

2.  Resulting  output  image  :  Fig.  B 
Step  2 

Apply  2-dimensional  3x3  low  pass  filter  to  original 
picture.  Fig.  B 

3.  Resulting  output  image  :  Fig.  C 
Step  3 

Perform  'AND’  operation  using  Fig.  B  as  input  'ONE'  to 
the  AND  operator  and  Fig.  C  as  input  'TWO'  to  the  AND 
operator . 

Threshold  =  0,  ie.  there  is  no  threshold  window 

4.  Resulting  output  image  :  Fig.  D 

Form  negative  image  of  Fig.  D  (ouput  image  from  step  3) 

5.  Resulting  ouput  image  :  Fig.  E 

Step  4 

Perform  'AND'  operation  using  Fig.  E  as  input  'ONE’  to 
the  AND  operator  and  Fig.  A  as  input  ’TWO'  to  the  AND 
operator . 

Case  1.  Threshold  =  (  the  Standard  Deviation  of  the 
pixel  values  for  the  entire  picture,  Fig.  A)  /  10 

6.  Resulting  ouput  image  :  Fig.  F 

Case  2.  Threshold  =  0,  ie.  no  threshold  window 

7.  Resulting  ouput  image  :  Fig.  E 

NOTE  :  Pixels  stored  in  output  image  of  AND  operation 
are  taken  from  input  ’TWO' 
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Other  images 


as  specified 


This  same  operation  is  repeated  exactly  for  a  1-Dimensional 
low  pass  filter.  The  only  change  in  the  processing 
procedure  occurs  in  Step  1  and  Step  2  as  follows* 

Step  1 

Apply  a  1-dimensional  1x7  low  pass  filter  to  original 
picture,  Fig.  A 

Step  2 

Apply  a  1-dimensional  1x3  low  pass  filter  to  original 
picture,  Fig.  A 

An  abbreviated  form  of  this  outline,  listing  only  the 
numbers  1  thru  8,  and  the  corresponding  correct  Figures,  is 
included  before  the  results  for  any  particular  image 
processed. 

Appendix  B  contains  a  three  dimensional  pixel  intensity 
(or  pi xel /pi cture  energy)  plot  for  each  of  the  original 
digitized  pictures  used  in  Appendix  A.  An  example  three 
dimensional  energy  plot  is  also  given  for  one  of  the 
Boolpass  Operator  1  output  images.  These  three  dimensional 
plots  help  demonstrate  and  emphasize  the  amount  of  picture 
energy  reduction  attained  through  use  of  a  Boolpass 
oper ator . 

During  the  research  and  testing  of  the  Boolpass 
operator  concept,  several  other  possible  useful  combinations 
for  Boolpass  operators  were  observed.  Time  and  resource 
limitations  did  not  allow  for  more  extensive  testing  of 
these  operators,  however,  the  initial  results  obtained  did 


-L. 
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warrant  -further  examination  -for  possible  future  development. 
Therefore,  Appendix  C  was  added  to  include  the  results  of 
several  of  these  other  possible  Boolpass  operators  and  may 
form  the  basis  for  future  research  and  testing. 

As  previously  noted,  the  examples  given  in  Appendix  A 
are  not  intended  to  represent  an  exhaustive  study  of  the 
Boolpass  operator  concept.  These  pictures  were  selected, 
and  their  Boolpass  operation  detailed,  because  (in  the  mind 
of  the  author)  they  appeared  to  have  the  greatest  potential 
for  use  in  an  image  processing/pattern  recognition  machine. 

For  clarity  and  ease  of  reference,  all  results  are 
presented  in  the  following  format: 

a)  Boolpass  operator  outline  (  of  image  processed) 

b)  digitized  image 

c)  intermediary  processed  images 

(in  sequential  order  of  processing) 

d)  final  output  image(s)  of  Boolpass  operator 

As  noted  previously,  pages  vi  thru  xviii  contain  a 
detailed  listing  of  all  pictures  displayed  in  Section  IV  - 
Results  and  Appendices  A  thru  C.  This  list  of  figures  may 
be  of  use  to  the  reader  as  an  additional  reference  aid. 
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V.  CONCLUSIONS  and  RECOMMENDATIONS 


Development  and  testing  of  Boolpass  operators,  as  a  new 
technique  for  real  image  processing,  was  presented  in  this 
thesis.  The  task  of  the  Boolpass  operators  was  one  of 
picture  energy  reduction,  while  retaining  the  fundamental 
picture  primitives  such  as  edges.  The  Boolpass  operators 
appeared  to  accomplish  this  energy  reduction  task,  with  a 
notable  degree  of  constancy,  over  a  number  of  very  different 
i mages. 

Five  factors  which  had  the  greatest  influence  during 
the  overall  development  and  testing  of  the  Boolpass 
operators  are  worthy  of  note.  These  factors  were  as 
foil owss 

1)  The  use  of  real  time  interactive  programming  techniques 
proved  to  be  a  powerful  tool  for  image  processing  operator 
development.  The  ability  to  see  (the  output  pictures), 
analyze,  and  respond  quickly,  throughout  all  stages  of  the 
processing,  was  instrumental  in  the  development  of  the 
Boolpass  operators. 

2)  The  nondestructive  processing  of  the  original  image  by 
use  of  additional  ouput  image  files  allowed  the  original 
image  to  remain  unchanged  for  further  processing  or 
reference  by  a  target  recognition  or  feature  extraction 
al gori thm. 

3)  In  those  pictures  where  multiple-levels  of  processing 
were  involved,  the  retention  of  the  image  greylevel  values, 
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as  opposed  to  going  to  a  binay  picture  representation, 
an  important  factor  in  the  'succesful'  application  of  a 
Boolpass  operator.  Also,  as  was  explained  in  Section  III, 
the  choice  of  which  input  image  pixel  values  should  be 
stored  in  the  output  image  file  can  have  a  profound  effect 
on  the  further  processing  of  that  image. 

4)  The  size  and  type  of  low  pass  filter  which  is  used  in 
processing  an  image  can  be  used  to  bias  an  image  for 
certain,  specific  features  or  orientations  of  those 
features.  This  was  illustrated  in  several  of  the  results 
where  only  one  dimensional  filtering  in  the  horizontal  axis 
was  used.  In  such  images,  the  Boolpass  operators  retained 
and  emphasized  horizonal  lines  while  supressing  many 
vertical  edges.  This  type  of  selective  filtering  could  have 
important  applications  for  aiding  target  recogntion. 
Similarly,  with  some  insight  into  a  particular 
pattern/target  recognition  problem,  design  or  choice  of  a 
low  pass  filter  specific  to  the  user's  requirements  could 
further  improve  the  Boolpass  operator  performance. 

5)  The  use  of  negative  images  is  an  important  step  which 
should  not  be  omitted  if  large  energy  reduction  is  desired 
through  use  of  the  Boolpass  operators  in  image  processing. 

This  thesis  has  presented  the  initial  development  and 
testing  of  a  new  technique,  the  Boolpass  operators,  for 
future  use  in  image  processing.  It  is  recommended  that 
further  testing  and  refinement  of  these  operators  be 
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continued  and  a  more  quantitative  evalutaion  of  their 
effectiveness  be  conducted.  More  specifically,  it  is 
recommended  that  further  research  include  at  least  these 
following  areas: 

1)  The  Boclpass  operators  should  be  exhaustively  tested  on 
a  very  wide  range  of  pictures,  with  large  variations  in 
pixel  mean  and  standard  deviation  for  the  pictures. 

2)  Other  types  of  thresholding  techniques  should  be 
investigated,  particularly  a  locally  adaptive  threshold 
techni que. 

3)  Different  types  of  low  pass  filtering  should  be 
exercised. 

4)  More  combinations  of  multilevel  processing  should  be 
attempted.  This  might  yield  useful  results  when  used  with  a 
different  filtering  technique. 

5)  The  performance  of  the  Bool  pass  operators  should  be 
quantitatively  compared  to  existing  edge  operators  such  as 
the  Sobel  and  Roberts  operators. 

6)  Research  into  the  possible  use  of  additional,  but 
different  operations  in  the  Boolpass  technique  should  be 
investigated.  For  example,  using  high  pass  filtering  of  the 
images  in  some  intermediary  stage  of  the  processing.  It 
might  also  be  found  that  other  Boolean  operations,  such  as 
exclusive  OR,  may  be  of  some  use  in  multiprocessing  the 
i mages. 

From  the  results  given  herein,  the  concept  and 
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techniques  of  the  Boolpass  operators  readily  demonstrate  a 
potential  for  future  development  and  possible  use  in  a 
pattern  recognition  machine. 
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Appendix  A 


Bool  pass  Ogeratgrs  _!  and  2  -  Results 


Contained  in  this  appendix  are  the  results  of  applying 
Boolpass  operators  1  and  2  to  eight  different  digitized 
images.  The  examples  which  follow  demonstrate  the  use  of 
both  one  and  two  dimensional  filtering  and  various  threshold 
values.  For  ease  of  reference,  the  general  outline  of  the 
processing  steps  for  both  Boolpass  operators  1  and  2  are 
repeated  herein. 

Outline  for  Bogl_gass  Operator  1 


1.  Digitized  Image  processed  :  Fig.  A 

Step  1 

Apply  2-dimensional  11x11  low  pass  filter  to  original 
picture,  Fig.  A 

2.  Resulting  output  image  :  Fig.  B 
Step  2 

Form  negative  image  of  Fig.  B  (ouput  image  from  step  1) 

3.  Resulting  ouput  image  :  Fig.  C 

Step  3 

Perform  ’AND’  operation  using  Fig.  C  as  input  ’ONE’  to 
the  AND  operator  and  Fig.  A  as  input  ’TWO’  to  the  AND 
operator 

Case  1.  Threshold  =  <1/2  the  Standard  Deviation  of 

the  pixel  values  for  the  entire  picture,  Fig.  A)  /  10 
A.  Resulting  ouput  image  :  Fig.  D 

Case  2.  Threshold  =  (the  Standard  Deviation  of  the 

pixel  values  fot  the  entire  picture,  Fig.  A)  /  10 
5.  Resulting  ouput  image  i  Fig.  E 
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NOTE  :  Pixels  stored  in  output  image  of  AND  operation 
are  taken  from  input  'TWO’ 


6.  Other  images  :  as  specified 


This  same  operation  is  repeated  exactly  for  a  1-Dimensional 
low  pass  filter.  The  only  change  in  the  processing 
procedure  occurs  in  Step  1,  as  follows; 


Step  1 

Apply  a  1-Dimensional  lxll  low  pass  filter  to  original 
picture,  Fig.  A 

An  abbreviated  form  of  this  outline,  listing  onlv  the 
numbers  1  thru  6  and  the  corresponding  correct  Figures  is 
included  before  the  results  for  any  particular  age 

processed . 


i 
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ABBendix  A  Qui^ck  Reference  to  BoDlgaes  QBgC®tQC5  i  §Qd  2 


Boolpass  Operator  1 


Page 

Fi gure 

2  dimensional  low  pass  filtering 

89 

43 

to 

to 

136 

83 

1  dimensional  low  pass  filtering 

137 

84 

to 

to 

160 

103 

Boolpass  Operator  2 


Page 

Figure 

2  dimensional  low  pass  filtering 

163 

104 

to 

to 

219 

153 

1  dimensional  low  pass  filtering 

220 

154 

to 

to 

267 

195 
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BOOLPASS  OPERATION  I 


1. 

Digitized  Image  Truck 

Fig.  43 

2. 

11x11  low  pass  filter 

Fig.  44 

3. 

Negative  of  2. 

Fig.  45 

4. 

AND  (1/2  Std  Dev) / 10 

Fig.  46 

5. 

AND  <  Std  Dev) / 10 

Fig.  47 

Mean : 

10  Standard  Deviation: 

4 
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BOOLPASS  OPERATION  I 


1. 

Digitized 

Image  Tomcat 

'• 

Fig.  48 

o 

X.  a 

11x11  low 

pass  -filter 

! 

Fig.  49 

3. 

Negati ve 

of  2. 

: 

Fig.  50 

4. 

AND  (1/2 

Std  Dev ) /  10 

a 

Fig.  51 

5. 

AND  (  Std 

Dev) / 10 

a 

a 

Fig.  52 

Mean:  10  Standard  Deviation:  4 
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I 


Fig.  51.  Output  Image  of  Boolpass  Operation  1 
Threshold  =  (1/2  Std  Dev)/10 
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BOOLPASS  OPERATION  I 


1 a 

Di  gi  t i s  ed 

Image  Field 

: 

Fig.  53 

o 

JL.  a 

11x11  low 

pass  -f i  1  ter 

• 

• 

Fig.  54 

3. 

Negati ve 

of  2. 

m 

Fig.  55 

4. 

AND  (1/2 

Std  Dev) / 10 

■ 

Fig.  56 

5. 

AND  (  Std 

Dev) / 10 

• 

a 

Fig.  57 

Means  8  Standard  Deviation:  5 
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Fi  g . 


54.  Field  Processed  by  an  Low  Pass  Filter 
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j; 


Fig. 


57.  Output  Image  of  Poolpass 
Threshold  =  <  Std  Dev  )/ 


Operation  1 
10 
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BOOLPASS  OPERATION  I 


1.  Digitized  Image  Field2  : 

2.  11x11  law  pass  filter  : 

3.  Negative  of  2.  : 

4.  AND  (1/2  Std  Dev > / 1 0  : 

5.  AND  <  Std  Dev) /10  : 

Mean:  9  Standard  Deviation:  5 


Fig.  61.  Output  Image  ot  Boolpass  Operation  1 
Threshold  =  <1/2  Std  Dev)/10 
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BOOLPASS  OPERATION  I 


1 . 

Digitized  linage  Hornet 

■ 

• 

Fig.  63 

2. 

11x11  low  pass  -filter 

■ 

Fig.  64 

3. 

Negative  of  2. 

Fig.  65 

4. 

AND  <1/2  Std  Dev) / 10 

Fig.  66 

5. 

AND  <  Std  Dev) / 10 

• 

■ 

Fig.  67 

Mean 

:  9  Standard  Deviation; 

3 
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Fi  g . 


67. 


Output  Image  of  Boolpass  Operation  1 
Threshold  =*  (  Std  Dev  )/  10 
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Fig, 


68. 


Digitised  Image  of  Original  Picture  Tank 


119 


121 


Fig.  71.  Output  Image  of  Boolpass  Operation  1 
Threshold  =  <1/2  Std  Dev)/10 
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BOOLPASS  OPERATION  I 

1.  Digitized  Image  Dispersal:  Pig.  73 

2.  11x11  low  pass  filter  :  Fig.  74 

3.  Negative  of  2.  :  Fig.  75 

4.  AND  <1/2  Std  Dev) / 10  :  Fig.  76 

5.  AND  <  Std  Dev) / 10  :  Fig.  77 

Mean:  9  Standard  Deviation:  5 


■er  i  isuvssan-  -rr^*r 

*V  il: 


/rsi 


-*./i 


X »“.*  i .'  . '  *  ”  ^  —  «l  ••J-.': ' 

V-;«'»£££  **,-  • 

'./•;  ~  :*  *  --nr- 

W%m^"  a 


J-~~:  *  ..•  ’  j  .-  ‘J  :"V" 

|  ••:•  ■•■  r.t:  •&/.•..,  ip-.-T-. 

£0 


•"■■■  !  M- 
> ::  *.-.•  •■■■■  •j-n* 


V-  ;••  i-5^;v*  r.  .s  “  %  ‘r-''i  '*•  *  \  HV 

■ ..  v  •"'  -A  -u  >  ^  * 

.  .-.■  i  .  -aa 

■■;••  ■■*  r.  g 

i  l  fv-“  *.  ''  ■'  "•■  ‘T  : 


itj 

:  ::£:::  r  i? 


:  &S 


’vWW 

Jsfesssilii 


73.  Digitized  Image  of  Original  Picture  Dispersal 


Fig. 


? 

i 
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74.  Dispersal  Processed  by  an  11x11  Low  Pass  Filter 
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BOOLF'ASS  OPERATION  I 


1 . 

Digitized  Image  Radar 

:  Fig.  78 

2. 

11x11  low  pass  -filter 

:  Fig.  79 

3. 

Negative  o-f  2. 

:  Fig.  80 

4. 

AND  (1/2  Std  Dev) / 10 

:  Fig.  81 

5. 

AND  (  Std  Dev) / 10 

:  Fig.  82 

6. 

AND  (Threshold  =  .4) 

:  Fig.  83 

Mean:  8 


Standard  Deviation:  5 
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BOOLFASS  OPERATION  I 


1 


1. 

Digitized  Image  Truck 

Fig.  84 

2. 

lxll  low  pass  filter 

■ 

« 

Fig.  85 

3. 

Negative  of  2. 

■ 

■ 

Fig.  86 

4. 

AND  <1/2  Std  Dev) /10 

■ 

■ 

Fig.  87 

5. 

AND  <  Std  Dev)/ 10 

» 

m 

Fig.  88 

Mean:  10  Standard  Deviation:  4 
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Fig.  87.  Output  Image  of  Boolpass  Operation  1 
Threshold  -  (1/2  Std  Dev)/10 
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Pig.  88.  Output  Image  of  Boolpass  Operation  1 
Threshold  =  <  Std  Dev  )/  10 
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BOOLPASS  OPERATION  I 


1. 

Digitized  Image  Tomcat 

! 

Fig.  99 

o 

A.  ■ 

lxll  low  pass  -filter 

• 

Fig.  90 

3. 

Negative  of  2. 

« 

• 

Fig.  91 

4. 

AND  (1/2  Std  Dev ) / 1 0 

• 

• 

Fig.  92 

5. 

AND  (  Std  Dev)/ 10 

■ 

■ 

Fig .  93 

Mean:  10  Standard  Deviation:  4 
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Fig.  89.  Digitized  Image  of  Original  Picture  F14  Tomcat 
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92.  Output  Image  of  Bool  pass  Operation  1 
Threshold  =  (1/2  Std  Dev>/10 
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Fig.  93.  Output  Image  of  Boolpass  Operation  1 
Threshold  =  <  Std  Dev  )/  10 
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BOOLPASS  OPERATION  I 

1.  Digitised  Image  Field2  : 

2.  lxll  low  pass  filter  : 

3.  Negative  of  2.  : 

4.  AND  <1/2  Std  Dev) / 10  : 

5.  AND  <  Std  Dev) / 10  : 

Mean:  9  Standard  Deviation:  E 
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g-  9 


ess  Operation 
Dev) / 10 


Output  Image  of  Boolpass  Operation  1 
Threshold  »  (  Std  Dev  )/  10 
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BOOLPASS  OPERATION  I 


1. 

Digitized  Image  Tank 

:  Fig.  99 

2. 

lxll  low  pass  filter 

:  Fig.  100 

3. 

Negative  of  2. 

:  Fig.  101 

A. 

AND  (1/2  Std  Dev ) / 1 0 

:  Fig.  102 

5. 

AND  (  Std  Dev) / 10 

:  Fig.  103 

Mean:  10  Standard  Deviation:  5 
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Fig.  99.  Digitized  Image  of  Original  Picture  Tank 
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102. 


Output  Image  o-f  Boolpass  Operati 
Threshold  =  <1/2  Std  Dev)/10 


15? 
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Outline  far  Bogl_gass  Oqeratgr  2 

1.  Digitized  picture  processed  :  Fig.  A 
Step  1 

Apply  2-dimensional  7x7  low  pass  filter  to  original 
pictur  Fig.  A 

2.  Resulting  output  image  :  Fig.  B 
Step  2 

Apply  2-dimensional  3x3  low  pass  filter  to  original 
picture,  Fig.  B 

3.  Resulting  output  image  :  Fig.  C 
Step  3 

Perform  'AND’  operation  using  Fig.  B  as  input  ’ONE’  to 
the  AND  operator  and  Fig.  C  as  input  ’TWO’  to  the  AND 
operator  Threshold  =  0,  ie.  there  is  no  threshold 
wi ndow 

4.  Resulting  output  image  :  Fig.  D 

Form  negative  image  of  Fig.  D  (ouput  image  from  step  3) 

5.  Resulting  ouput  image  :  Fig.  E 

Step  4 

Perform  ’AND'  operation  using  Fig.  E  as  input  ’ONE’  to 
the  AND  operator  and  Fig.  A  as  input  ’TWO’  to  the  AND 
operator 

Case  1.  Threshold  =  <  the  Standard  Deviation  of  the 

pixel  values  for  the  entire  picture.  Fig.  A)  /  10 

6.  Resulting  ouput  image  :  Fig.  F 

Case  2.  Threshold  =  0,  ie.  no  threshold  window 

7.  Resulting  ouput  image  :  Fig.  E 

NOTE  :  Pixels  stored  in  output  image  of  AND  operation 
are  taken  from  input  ’TWO' 

8.  Other  images  :  as  specified 

This  same  operation  is  repeated  exactly  for  a  1-Dimensional 
low  pass  filter.  The  only  change  in  the  processing 
procedure  occurs  in  Step  1  and  Step  2  as  follows: 


161 


Step  1 

Apply  a  1-Dimensional  1x7  low  pass  filter  to  original 
picture,  Fig.  A 

Step  2 

Apply  a  1-dimensional  1x3  low  pass  filter  to  original 
picture,  Fig.  A 

An  abbreviated  form  of  this  outline,  listing  only  the 
numbers  1  thru  8,  and  the  corresponding  correct  Figures,  is 
included  before  the  results  for  any  particular  image 
processed . 
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BOOLPASS  OPERATION  2 


1. 

Digitized  image  Truck 

:  Fig. 

104 

2. 

7x7 

low  pass  filter 

:  Fig. 

105 

3. 

3x3 

low  pass  filter 

:  Fig. 

106 

4. 

AND 

(0  Std  Dev) 

:  Fig. 

107 

5. 

Negative  of  4. 

:  Fig. 

108 

6. 

AND 

(Std  Dev ) / 1 0 

:  Fig. 

109 

7. 

AND 

(0  Std  Dev)/10 

:  Fig. 

110 

Mean:  10 


Standard  Deviation:  4 


107. 


Boolean  AND  Operation  Performed  On  Prereeding 
Two  Figures.  Threshold  =  0 


167 


•  i.  r  •  •  r  r*.  t  ’  i-r . -w.**  r  )  ..  >rj>  msmjn  -fvrvjv,;.  =  .u-.  .  •  — , .  f ,  ,  rv,  •*>  ;;  v.r  - ;:  -4 

**•  f.  •  **«  »** ’*  t.;*' '?K  i5  #  'f-  •  ;’■*'*  ****,IS  j-2.|r»  i  :*  *  -j  •“  f  *-•  i  ^  :-  J  r*  ‘  ^  .  ; 

r>-:\  I'y  i»r. ;-h  rv;  'irig ^.^iS «;•:>•■  v?‘ f.!  ■ >.■!! i!!v  •! lV •!-?£ 

.  4—  *  f-J-iiltef  »sil  I&  1  ^  ?.  ?J  ’•  •;  r  ii- :::  ::*3l 

K  -  S.'-~  ?#J. f::v  *r*  y*.-i  ’*??  *  J#  ? •*  -  k  •  v,**;*-^ 


«.  1^; 


•Jr--*  .  •• 

v*  ■  - 


■'»  p; ;-  r~  -  i  s,  "i r  J? ; .“  \  v»  ':  •'  -f  >»  il  rsH  ■ .-  >i  .i  :•. k;  v  V:-  < 

K,f  *  -  "  •  -ri  •  • ♦  '  *' .!'.*•  •  *  *|  T»*-«  ,1  **i  -  't  **>-.•  .r.  •  ^  ••'..»  i-  ;j;  }  f  *  tf  9.  J  f  -2* 

'-•j-r  ;'}■>  ir  A..'-  ,  If  ■■«•>■  .V|  ^..:J  -  .••  .:  •..  ui^i  t  -i  'M 

'  •%'  ‘  V*;  •VT.*'  t’  'i-W  ’ _  —  ”/ "#  f-  i“  rV.  •  •  V' I .tit.  ri:  • 


Fig.  108.  Negative  Image  of  Fig.  107 
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Output  Image  O-f  Bool  pass  Operation  2 
Threshold  =  <  Std  Dev  )/10 


BOOLPASS  OPERATION 


1 . 

Digitized  image  Tomcat 

Fig. 

1 1 1 

o 

X-  m 

7x7 

low  pass  filter 

Fig. 

1 12 

3. 

3x3 

low  pass  filter 

■ 

• 

Fig. 

113 

4. 

AND 

(0  Std  Dev) 

■ 

Fig. 

114 

5. 

Negative  of  4. 

Fig . 

1  15 

6. 

AND 

( Std  Dev ) / 1 0 

Fig. 

116 

7. 

AND 

<0  Std  Dev) / 1 0 

Fig. 

117 

Mean 

:  10 

Standard  Deviation: 

4 
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Fig.  114.  Boolean  AND  Operation  Performed  On  Preceeding 
Two  Figures.  Threshold  =  0 
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115.  Negative  Image  of  Fig.  114 
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BOOLPASS  OPERATION  2 


1. 

Digitized  image  Field 

s 

Fig. 

118 

■£_  • 

7x7  low  pass  filter 

: 

Fig. 

119 

3. 

3x3  low  pass  filter 

■ 

■ 

Fig. 

120 

■4. 

AND  <0  Std  Dev) 

• 

m 

Fig. 

121 

5. 

Negative  of  4. 

m 

■ 

Fig. 

122 

6. 

AND  <Std  Dev) / 10 

0 

* 

Fig. 

123 

7. 

AND  <0  Std  Dev) /10 

: 

Fig. 

124 

Mean:  Q  Standard  Deviation:  5 
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120.  Field  Processed  BY  A  3x3  Low  Pass  Filter 


1 82 


Fig.  122.  Negative  Image  of  Fig.  121 
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Fig.  123.  Output  Image  Q-f  Boolpass  Operation  2 
Threshold  =*  <  Std  Dev  )/10 
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BOOLPASS  OPERATION  2 


1. 

Digitized  image  Field2 

Fig. 

125 

2. 

7x7  low  pass  -filter 

Fig. 

126 

3. 

3x3  low  pass  -filter 

Fig. 

127 

4. 

AND  (0  Std  Dev) 

• 

• 

Fig. 

128 

5. 

Negative  of  4. 

• 

■ 

Fig. 

129 

6. 

AND  <  Std  Dev) / 1 0 

■ 

■ 

Fig. 

130 

7. 

AND  (0  Std  Dev) /10 

s 

Fig. 

131 

Mean:  9  Standard  Deviation:  5 
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Fig.  12S.  Boolean  AND  Operation  Performed 
Two  Figures.  Threshold  *  0 


On  Pr eceeding 


191 
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Fig.  131.  Output  Image  Of  Boolpass  Operation  2 
Threshold  =  0 
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BOOLPASS  OPERATION 

2 

1 . 

Digitized  image  Harriet 

• 

• 

Fig.  1 32 

2. 

7x7  low  pass  -filter 

m 

• 

Fig.  133 

3. 

3x3  low  pass  filter 

■ 

■ 

Fig.  134 

4. 

AND  (0  Std  Dev) 

• 

■ 

Fig.  135 

5. 

Negative  of  4. 

: 

Fig.  136 

6. 

AND  (Std  Dev) 710 

■ 

■ 

Fig.  137 

7. 

AND  (0  Std  Dev) 710 

: 

Fig.  138 

Standard  Deviation:  3 
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Fig.  133.  F18  Hornet  Processed  By  A  7x7  Low  Pass  Filter 
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Boolean  AND  Operation  Performed  On  Preceeding 
Two  Figures.  Threshold  =  0 
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Fig.  137.  Output  Image  Of  Boolpass  Operation  2 
Threshold  =  (  Std  Dev  )/10 
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BOOLPASS  OPERATION  2 


1. 

Digitized  image  Tank 

m 

• 

Fig. 

139 

2. 

7x7  low  pass  filter 

m 

• 

Fig. 

140 

3. 

3x3  low  pass  filter 

• 

• 

Fig. 

141 

4. 

AND  (0  Std  Dev) 

• 

* 

Fig. 

142 

5 . 

Negative  of  4. 

: 

Fig. 

143 

6. 

AND  (Std  Dev) / 10 

■ 

Fig. 

144 

7. 

AND  (0  Std  Dev) / 10 

• 

• 

Fig. 

145 

8. 

AND  (Threshold  =  .4) 

• 

• 

Fig. 

146 

Mean:  10  Standard  Deviation:  5 
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205 
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BOOLPASS  OPERATION  2 

1.  Digitized  image  Dispersals  Fig.  147 

2.  7>;7  low  pass  -filter  :  Fig.  148 

3.  3x3  low  pass  filter  :  Fig.  149 

4.  AND  <0  Std  Dev)  :  Fig.  150 

5.  Negative  of  4.  :  Fig.  151 

6.  AND  < Std  Dev) / 10  :  Fig.  152 

7.  AND  (0  Std  Dev)/10  s  Fig.  153 

Mean:  9  Standard  Deviation:  5 
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Fig 


.  153. 

I 


Output  Image  Of  Boolpass  Operation  2 
Threshold  *  0 


+  ttt 


Digitized  video  image  of  Truck 


Truck 


Fi 1  ter 


Fig.  156.  Truck  Processed  BY  A  1x3  Low  Pass  Filter 
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224 


Fig. 


159.  Output  Image  Of  Bool  pass  Operation 
Threshold  =  (  Std  Dev  >/10 
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Output  Image  Of  Boolpsss  Operat 
Threshold  =  0 


BOOLPASS  OPERATION  2 


Digitized  image  Tomcat  :  Fi 
1x7  low  pass  -filter  s 

1>:3  low  pass  filter  s 

AND  <0  Std  Dev)  s 

Negative  of  4.  : 

AND  (Std  Dev) / 10  : 

AND  <0  Std  Dev ) / 1 0  : 

10  Standard  Deviation 


A 


Fig.  162.  Tomcat  Processed  By  A  1x7  Low  Pass  Filter 
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I 

l-_ 
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Fig . 


163.  Tomcat  Processed  BY  A  1x3  Low  Pass  Filter 
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Fig.  164.  Boolean  AND  Operation  Performed  On  Preceeding 
Two  Figures.  Threshold  =  0 
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165.  Negative  Image  of  Fig. 


233 


167. 


Output  Image  Of  Boolpass  Operation  2 
Threshold  ®  0 
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BOOLPASS  OPERATION  2 


1. 

Digitized  image  Field 

5 

Fig. 

16B 

2. 

1x7 

low  pass  -filter 

Fig. 

169 

3. 

1x3 

low  pass  filter 

* 

Fi  g. 

170 

4. 

AND 

<0  Std  Dev) 

• 

Fig. 

171 

5. 

Negative  of  4. 

: 

Fig. 

172 

6. 

AND 

<Std  Dev) /10 

■ 

• 

Fig. 

173 

7. 

AND 

<0  Std  Dev) /10 

• 

Fig. 

174 

Mean:  8  Standard  Deviation:  5 
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172.  Negative  Image  of  Fig.  171 


i' 


Fig.  173.  Output  Image  Of  Boolpass  Operation 
Threshold  =  <  Std  Dev  )/10 
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BOOLPASS  OPERATION  2 


Digitized  image  Field2 

■ 

Fig.  3 

1x7  low  pass  filter 

■ 

• 

Fig.  1 

1x3  low  pass  filter 

■ 

■ 

Fig. 

AND  (0  Std  Dev) 

a 

Fig. 

Negative  of  4. 

a 

• 

Fig. 

AND  <Std  Dev) / 10 

a 

■ 

Fig. 

AND  (0  Std  Dev)/ 10 

a 

Fig. 

9 


Standard  Deviation:  5 
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Fig.  179.  Negative  Image  of  Fig.  178 
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Fig 


181 


Output  Image  Of  Bool  pass  Operation  2 
Threshold  =  0 
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BOOLF'ASS  OPERATION  2 


1 . 

Digitized  image  Tank 

Fi  g. 

182 

2. 

1x7  low  pass  filter 

• 

• 

Fig. 

183 

3. 

1x3  low  pass  filter 

• 

Fig. 
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4. 

AND  (0  Std  Dev) 

Fig. 
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5. 

Negative  of  4. 

• 

• 

Fig. 
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6. 

AND  (Std  Dev) /10 

• 

• 

Fig. 

187 

7. 

AND  (0  Std  Dev) / 10 

• 

■ 

Fig. 

188 

Mean 

:  10  Standard  Deviation: 

5 

Digitized  video  image  of  Tank 


Fig 


183.  Tank  Processed  By  A  1x7  Low  Pass  Filter 
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Fig.  187.  Output  Image  Of  Boolpass  Operation  2 
Threshold  =  (  Std  Dev  )/10 
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BQOLPASS  OPERATION  2 


1. 

Digitized  image  Dispersal 

s 

Fig. 

189 

2. 

1x7  low  pass  filter 

Fig. 

190 

3. 

1x3  low  pass  filter 

• 

a 

Fig. 

191 

4. 

AND  (0  Std  Dev) 

• 

• 

Fig. 

192 

5. 

Negative  of  4. 

• 

a 

Fig. 

193 

6. 

AND  (Std  Dev) / 10 

• 

■ 

Fig. 

194 

7. 

AND  (0  Std  Dev) /10 

a 

a 

Fig. 

195 

Mean:  9  Standard  Deviation:  5 
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191.  Dispersal  Processed  BY  A  1x3  Low  Pass  Filter 
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192.  Boolean  AND  Operation  Performed  On  Preceeding 
Two  Figures.  Threshold  =  0 
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Fig 


194.  Output  Image  Of  Boolpass  Operation  2 
Threshold  *  (  Std  Dev  >/10 

i 
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COMBINED  SPATIAL  FILTERING  AND  BOOLEAN  OPERATORS  APPLIED  TO  THE — ETCUI) 
JUN  82  B  E  FELTMATE 

UNCLASSIFIED  AFIT/GCS/EE/82J-8  NL 


AgeeQdix  B 

SnQLCtDSLQQlL  Pi  cture/Pi.xel_  Intensity  Plots 

The  following  3-dimensional  plots  help  to  visually 
display  to  the  reader  the  exact  pixel  intensity  distribution 
over  an  entire  digitized  image.  This  appendix  contains 
pixel  intensity  plots  for  each  of  the  original  digitized 
images  used  in  Appendix  A.  An  example  3-dimensional 
intensity  plot  is  also  given  for  one  of  the  Boolpass 
Operator  1  output  images  (for  each  of  the  eight  original 
images  used  in  Appendix  A). 

A  large  amount  of  overall  energy  reduction  is  apparent 
in  all  images  processed  by  Boolpass  Operator  1. 

NOTE:  Anomalies  existed  in  the  computati onal  routines  used 

for  making  these  plots.  Therefore,  the  graphs  which  follow 
are  reversed  from  the  original  image.  For  example,  an 
aircraft  pointed  toward  the  right  side  of  the  original  image 
points  to  the  left  side  of  the  image  after  being  processing 
by  the  3-dimensional  plot  routines.  Despite  this 
distraction,  the  plots  are  still  considered  constructive 
visual  aids  to  the  better  understanding  of  the  image 
processing  operation  done  by  a  Boolpass  operator. 
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t  mum=  12.0000 


Fig. 


199. 


3-Dimensional  Pixel  Plot  of  Fig. 
Bool  pass  Operator  1  Used  for  Pro 
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Fig. 


201. 


3-Dimensional  Pixel  Plot  of  Fig. 
Bool  pass  Operator  1  Used  For  Pr oi 
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Fig.  202 


3-Di mensi onal 


3-Dimensional  Pixel  Plot  of  Fi 
Boolpass  Operator  1  Used  For 


Fig.  208.  3-Dimensional  Pixel  Plot  of  D 
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Additional  Boglgass  Ogeratgrs  —  Results 

The  number  of  combinations  of  low  pass  filtering 
techniques  and  boolean  operations  which  might  constitute  a 
Boolpass  operator  are  very  large.  Appendix  A  demonstrated 
two  specific  combinations  which  were  labelled  Boolpass 
operators  1  and  2.  These  operators  were  then  applied  in  a 
consistent  manner  to  8  different  digitized  images  so  that 
relative  comparisons  of  the  output  images  were  possible. 

The  results  given  in  this  section  demonstrate  other 
possible  Boolpass  operators  that  were  briefly  investigated 
during  the  research  and  development  of  this  thesis.  These 
additional  Boolpass  operators  were  chosen  for  inclusion  in 
this  appendix  based  on  their  potential  for  possible  future 
testing  and  refinement.  At  this  time,  however,  they  have 
not  been  extensively  tested  over  a  large  number  of  images. 

Figures  212  thru  216  are  variations  on  Boolpass 
Operation  1,  as  previously  outlined  in  Appendix  A.  The 
differences  in  the  processing  steps  for  these  six  images 
are  detailed  below.  Only  the  output  images  from  the 
resulting  Boolpass  operation  are  given  in  this  section. 
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y*Ci stiQQs  On  Bool_gass  Operation  1 

Figure  212: 

Original  digitized  Image  -  Field  <Fig.  53).  Image 
Field  is  replaced  (as  input  TWO  to  the  AND  operation)  by  3x3 

low  pass  -filtered  image  of  Field  (Fig.  120). 

Threshold  =  .2 

Figure  213: 

Original  digitized  image  -  Field  (Fig.  53).  Image 
Field  is  replaced  (as  input  Two  to  the  AND  operation)  by  7x7 

low  pass  filtered  image  of  Field  (Fig.  119). 

Threshold  =  .2 

Figure  214: 

Original  digitized  image  -  Field  (Fig.  53).  11x11  low 

pass  filtered  image  of  Field  is  replaced  (as  input  ONE  to 
the  AND  operation)  by  a  7x7  low  pass  filtered  image  of  Field 
(Fig.  119) 

Threshold  =  .5 

Figure  215: 

Original  digitized  image  -  Field2  (Fig.  58) .  11x11  low 

pass  filtered  image  of  Field2  is  replaced  (as  input  ONE  to 
the  AND  operation)  by  a  7x7  low  pass  filtered  image  of 
Field2  (Fig.  126).  Similarly,  image  Field2  is  replaced  (as 
input  Two  to  the  AND  operation)  by  5x5  low  pass  filtered 
image  of  Field2. 

Threshold  =  .3 
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Figure  216: 

Original  digitized  image  Dispersal  (Fig.  73).  Image 
Dispersal  is  replaced  (as  input  Two  to  the  AND  operation)  by 
7x7  low  pass  -filtered  image  o-f  Dispersal  (Fig.  14B)  . 

Threshold  =  .2 

Figure  217: 

Original  digitized  image  Truck  (Fig.  7).  Image  Truck 
is  replaced  (as  input  Two  to  the  AND  operation)  by  7x7  low 
pass  filtered  image  of  Truck. 

Threshold  =  .4 

Figure  218: 

Original  digitzed  image  Tomcat  (Fig.  48).  Image 
Tomcat  is  replaced  (  as  input  Two  to  the  AND  operation)  by 
7x7  low  pass  filtered  image  of  Tomcat. 

Threshold  =  0 

Figure  219:  Threshold  =  .2 

Figure  220:  Threshold  =  .4 

Figure  221: 

Original  digitized  image  Dispersal  (Fig.  73)  1x5  low 

pass  filtered  image  of  Dispersal  is  AND'ed  with  a  1x7  low 
pass  filtered  image  of  Dispersal.  Threshold  for  this  first 
AND  operation  =  .3  .  The  ouput  image  of  this  first  AND 
operation  is  used  as  input  ONE  to  a  second  AND  operation, 
and  Dispersal  (the  original  digitized  image)  is  used  as 
input  Two  to  the  AND  operation.  Figure  221  is  the  ouput  of 
this  second  and  operation  performed  at  a  threshold  *=  . 
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Combined  Images  Using  Boolean  OR  Operation 


Figures  222  and  223  are  each  composed  of  a  combination 
of  two  output  images  from  Boolpass  operations.  The 
combining  of  the  two  images  was  performed  by  a  simple 
boolean  OR  operation.  For  the  two  examples  which  follow, 
one  image  for  the  OR  operation  was  an  ouput  image  from 
Boopass  operation  1,  while  the  second  image  for  the  OR 
operation  was  an  ouput  image  from  Boolpass  operation  2. 

Figure  222: 

Original  digitized  image  -  Field  (Fig.  53).  Boolpass 
operation  1  output  image,  Fig.  213  OR’ed  with  Boolpass 
operation  2  ouput  image.  Fig.  124. 

Figure  223: 

Original  digitized  image  -  Dispersal  (Fig.  73). 
Boolpass  operation  1  output  image,  Fig.  76  OR’ed  with 
Boolpass  operation  2  ouput  image.  Fig.  153. 

Bgolpass  Operation  3 

Because  Figure  224  uses  a  more  complex  procedure  for 
processing,  a  complete  outline  of  the  processing  steps  is 
included  for  this  image. 
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Outline  for  Bgolpass  Operator  3 


1.  Digitized  picture  processed  :  Fig.  A 
Step  1 

Apply  2-dimensional  3x3  low  puss  filter  to  original 
picture,  Fig.  A 

2.  Resulting  output  image  :  Fig.  B 
Step  2 

Apply  2-dimensional  7x7  low  pass  filter  to  original 
picture,  Fig.  A 

3.  Resulting  output  image  :  Fig.  C 
Step  3 

Perform  ’AND’  operation  using  Fig.  B  as  input  ’ONE’  to 
the  AND  operator  and  Fig.  C  as  input  ’TWO’  to  the  AND 
operator  Threshold  =  0,  ie.  there  is  no  threshold 
wi ndow 

4.  Resulting  output  image  :  Fig.  D 
Step  4 

Apply  2-dimensional  5x5  low  pass  filter  to  original 
picture.  Fig.  A 

5.  Resulting  output  image  :  Fig.  E 
Step  5 

Form  negative  image  of  Fig.  D  (ouput  image  from  step  3) 

6.  Resulting  ouput  image  :  Fig.  F 

Step  6 

Perform  ’AND’  operation  using  Fig.  E  as  input  'ONE'  to 
the  AND  operator  and  Fig.  F  as  input  'TWO'  to  the  AND 
oper ator 

Threshold  =  (  the  Standard  Deviation  of  the  pixel 
values  for  the  entire  picture.  Fig.  A)  /  10 

7.  Resulting  ouput  im^ge  :  Fig.  6 

NOTE  :  Pixels  stared  in  output  image  of  AND  operation 
are  taken  from  input  ’TWO' 
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Variation  on  Boolpass  Operation  1.  Image  Field 
is  Replaced  (as  input  TWO  to  the  AND  operation) 
by  3x3  Low  Pass  Filtered  Image  of  Field  (Fig. 
120).  Threshold  =  .2 


Fig.  213.  Variation  on  Boolpass  Operation  1.  Image  Field 
is  Replaced  (as  input  Two  to  the  AND  operation) 
by  7x7  Low  Pass  Filtered  Image  o-f  Field  (Fig. 
119).  Threshold  =■  .2 
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Variation  On  Bool  pass  Operation  1.  11x11  Low 
Pass  Filtered  Image  of  Field  is  Replaced  (as 
input  ONE  to  the  AND  operation)  by  a  7x7  LPF 
Image  of  Field  (Fig.  119).  Threshold  ■  .5 


Fig.  215.  Variation  On  Boolpass  Operation  1.  11x11  LPF 
Image  of  Field2  is  Replaced  (as  input  ONE  to  the  AND 
operation)  By  a  7x7  LPF  Image  of  Field2  (Fig.  126).  Image 
Field2  is  Replaced  (as  input  Two  to  the  AND  operation)  by 
5x5  LPF  Image  of  Field2.  Threshold  =  .3 


Fig.  217.  Original  digitized  image  Truck  (Fig.  7).  Image 
Truck  is  replaced  (as  input  Two  to  the  AND 
operation)  by  7x7  low  pass  filtered  image  of 
Truck.  Threshold  =  .4 
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Fig.  218.  Original  digitred  image  Tomcat  <Fig.  48).  Image 
Tomcat  is  replaced  <  as  input  Two  to  the  AND  operation)  by 
7x7  low  pass  filtered  image  of  Tomcat.  Threshold  “  0 
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Fig.  222.  Original  Digitized  Image  -  Field  (Fig.  53). 

Boolpass  Operation  1  Output  Image,  Fig.  213 
OR’ed  With  Boolpass  Operation  2  Output  Image, 
Fig.  124. 
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Fig. 


OR’ed  With  Boolpass  Operation  2  Output  Image, 
Fig.  153. 
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Appendix  D 


CQCByter  Program  Listings 

The  computer  programs  used  for  image  processing  with 
the  Boolpass  operators  performed  two  main  functions:  1)  low 
pass  filtering,  and  2)  boolean  operations.  The  programs, 
as  presently  written,  are  general  in  nature  and  have  not 
been  performance  optimized  to  any  Boolpass  operator  in 
particular.  It  should  also  be  noted  that  this  RDQS  version 
of  Pascal  was  limited  to  a  maximum  stack  size  of  16K, 
imposing  large  disk  I/O  times  for  the  processing  of  64K 
digitized  images. 

In  order  to  overcome  the  limitations  imposed  by  the 
maximum  stack  size,  the  64K  video  disk  files  were  divided 
into  four  equal  files  of  1 6K  each  for  subsequent  processing. 
After  processing,  these  four  files  were  then  sequentially 
appended  back  into  a  64K  video  file.  This  latter  step  was 
required  to  be  able  to  display  the  picture  file  on  the  AFIT 
signal  processing  laboratory  video  monitor. 

Two  Fortran  programs  were  used  during  the  processing  of 
the  Boolpass  operators.  One  program  was  used  for  dividing 
the  64K  video  files  into  four  16K  files.  The  second  Fortran 
program  was  used  for  creating  the  negative  of  a  video  file. 
These  programs  already  existed  as  utility  programs  in  the 
AFIT  signal  processing  1 aboratory  and  are  not  included  in 
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,  *  MIMM  «««  ««•»»»»**»*  «***»*#»*  D  »  «  M  •«  «  «  »  «  « 

PROGRAM  NAME  AVGMA5KS  PA 

WRITTEN  BY  CAPTAIN  BLAINE  FElTMATE 

DATE  SEP  1981  REVISED  21  MAY  1982 


PURPOSE  the  PURPOSE  OF  THIS  PROGRAM,  AVCMASKS  PA,  IS  TO  ACT  AS  A  LOUPASS 

FILTER  IN  PROCESSING  DIGITIZED  VIDEO  IMAGES  THE  FILTER  IS  A  SIMPLE  AVERAGING 
OF  THE  MASK  PIXEL  VALUES,  WITH  THE  CENTER  PIXEL  BEING  REPLACED  BY 
THE  AVERAGE'  VALUE  .  THIS  PROGRAM  CAN  WORK  WITH  A  MAXIMUM  PICTURE  SIZE  OF 
256  X  256  PIXEL  ELEMENTS,  AND  A  MAXIMUM  OF  16  CREYLEVEL  VALUES. 

DUE  TO  THE  STACK  LIMITATIONS  IMPOSED  BY  THE  RDOS  VERSION  OF 
PASCAL.  THE  64K  PICTURE  FILE  MUST  BE  EQUALLY  DIVIDED  INTO  FOUR  SEPARATE 
FILES  NAMED  DATA1, DATA2, DATA3. DATA4  AFTER  PROCESSING  BY  THIS  PROGRAM 
The  new  PICTURE  WILL  BE  STORED  IN  FILES  NEWD1,  NEWD2,  NEWD3,  NEUD4 
IN  ORDER  TO  DISPLAY  THE  PICIvRE  ON  THE  MONITOR  VIA  THE  SIGNAL  PROCESSING 
LAB  PROGRAM  "VIDEO”  THE  SEPARATE  FILES  MUST  AGAIN  BE  RECOMBINED  INTO  ONE 
SINGLE  PICTURE  FILE.  THIS  IS  EASILY  DONE  BY  THE  FOLLOWING  COMMAND 
"  APPEND  <FILENAME>  NEWD1  NEWD2  NEWD3  NEWD4" 

TO  INITIALLY  DIVIDE  THE  PICTURE  FILE  INTO  FOUR  -  USE  THE  FORTRAN  PROGRAM 
LABELLED  “V4M.  BY  RUNNING  V4  YOU  WILL  BE  PROMPTED  FDR  THE  PICTURE 
FILENAME,  AND  ONCE  THE  FILENAME  IS  ENTERED  EVERYTHING  ELSE  PROCEEDS  AUTOMAT¬ 
ICALLY  WITHIN  SECONDS 

AT  THIS  TIME  THE  SIZE  MASKS  ALLOWED  FOR  THIS  PROGRAM  ARE  3X3.  5X5,  7X7. 

9X9.  11X11 

IF  LARGER  MASKS  SI7ES  ARE  REQUIRED  MODIFICATIONS  TO  THE  PROGRAM  WILL 
BE  REQUIRED  -  (  BUT  COULD  BE  VERY  EASILY  DONE). 


am  «*««*  ***•  #  ##«*««  #»««#«  »«#«««  #  ««*««  #»****##«»  *«*  *«»#«  ««««*  ««**#«  •  > 


PROGRAM  AVGMASKSt INPUT.  OUTPUT.  DAT  A  1 ,  DAT  A2> DAT A3.  DAT  A4, 
NEWD1.  NFWD2.  NEWD3. NEWD4); 


LABEL  e2. 


YfSCHAR  =  'Y'i 
NOCHAR  «  'N'. 


T  T  PE 

DATAFILE  =  PACKED  RECORD 

LINE  .  PACKED  ARRAY! 0  255]  OF  0  15 

END; 

DAT AF  =  FILE  OF  DATAF II E» 

ROWCOL  *  O.  .  255. 

PIXEL  -  0  15; 

RfcCROWDATA  *  PACKED  ARRAY! 1 .  11]  OF  DATAFILE; 

RECP1XM  »  PACKED  RECORD 

PIXMATRIX  .  PACKED  ARRAY  C 1 .  121]  OF  0  15; 

END, 

RECMATRI  XMA5K  =*  PACKED  ARRAY  [  0  255]  OF  RECF  I  XM. 


VAR 

GRECMASK  RECMATR I XMASK; 

MATRnrDATA  RECROWDATA. 

DATAOBJECT  DATAFILE. 

DATA1 , DATA2.  DATA3- DATA4  DATAF,  ( * THESE  ARE  THE  INPUT  IMAGE  FILES  TO  BE  PROCESSED*) 

NEWD1 . NEWD2.  NEWD3.  NEWD4  DATAF,  < *THESE  ARE  THE  OUPUT  IMAGE  FILES  WHICH  CONTAIN  THE  PROCESSED  IMAGE*) 
SREYLEVEL  PIXEL, 
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INTEGER. 


( • 


=iOw.  COu  ROwCQL* 
f-AXCOL.  MjNCOL-  MAXROW.  HINROW 
ANSWER  ChaR, 

LINECOUNT,  *oun 7.  mask,  number  integer. 

OFFSET. SI2E  INTEGER, 

MEAN. STDEVIAT ION  INTEGER. 


*l***##«l»****<t#*> 


..PROCEDURE  PLANwFILE  IS  DESIGNED  TO  PILL  THE  OUTPUT  DISK 
NEWD1  AND  NEWD2  WITH  THE  LOWEST  GREYLEVEL  VALUE. 1  E 
■15'  WHICH  IS  EQU I  VALENT  TO  A  BLANK  THIS  tNSURES  .HAT  NO 
IS  CONTAINED  IN  THESE  FILES  BEFORE  PICTURE  ELEHENT  INFORHAT 


FILES 

'GARBAGE ' 

ION  IS  STORED*) 


PROCEDURE  BLANKFILE. 


VAR 


num.  index 

N1PICtURE 

N2P1CTURE 


0  255. 

packed  ARRAVCO  633  OF  DATAFILE. 
PACKED  ARRAYC64  127]  QF  DATAFILE. 


RFWR I TE (Nr WD1 ) . 

FOR  INDEX  *  0  TO  63  DO 
BEGIN 

WITH  NlPIOUREt  INDEX]  DO 
BEGIN 

FOR  NUM  '0*0  255  DO 
LINE! NUM]  =  15. 

END. 

NE  WD 1  **  =  NlPlCTl'RE  [  INDEX  3* 

PUT  <  NEWD 1 ) - 


END. 

Cl  USE  <  NEWDJ : . 

RE WR 1 T  E  <  NFWD2 ) . 

FOR  INDEX  *  64  TO  1 
BEGIN 

WITH  N2PICTUREC 
BEGIN 

F  OR  NUM 
L 1NEI NUM] 


27  DO 

INDEX]  DO 

*  0  TO  255  DO 
=  15. 


END. 

NF WD2 **  =N2P  ICTUREC  INDEX]. 

PUT (NFWD2) . 


END. 

CLOSE ( NFWD2 ) i 

END. 


(  .BLANK  OF  ILE  IS  DESIGN.  3  TO  FILL 
WITH  THE  LOWEST  GREYlEVEL  VALUE. 


THE  OUTPUT  DISK  FILES  NEWD3  AND 
I  E  '15'.  WHICH  IS  EQUIVALENT 


NEWD4 

TO  A  BLANK 


•  ) 


PROCEDURE  BLANK2FILE. 


VAR 

Nun. INDEX  0  ?55. 
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N3PJC71/RE  PACKED  ARRAY112B  19]  3  OP  DATAFILE. 
N4P2CTURE  PACKED  ARRA>[192  2553  OF  DATAFILE, 

BEGIN 

REWRI TE (NEWD3) , 

FOR  INDEX  .=  128  TO  191  DO 
BEGIN 

WITH  N3P ICTUREC INDEX!  DO 
BEGIN 

FOR  NUM  =  C  TO  255  DO 
LINE  [  NUM  3  =  15, 

END. 


END 


NEUD3~  .  *  N3P1CTUREC INDEX  3. 
PUT  *  NEWD3) ; 

END, 

CLOSE <  NEWD3 ) ; 

REWR I TE ( NE WD4 )  ; 

FOR  INDEX  : =  192  TO  255  DO 
BEGIN 

WITH  N4P ICTUREC INDEX!  DO 
BEGIN 

FOR  NUM  =  0  TO  255  DO 
LINE  C NUM!  .  *  15, 

END. 

NEb'D4A  *■  N4P  ICTUREE  INDEX  3. 
PUT ( NL WD4 )  i 

END; 

CLOSE (NEWD4). 

£ •  PROCEDURE  BLANKFILEO 


f  »  tM>  «»»«««  »  S 


<  * PROCEDURE  Cm^CKLIMITS  IS  DESIGNED  TO  CHECK  THE  ROW  AND  COLUMN  LIMITS 
WHICH  ARE  USED  f-OR  PROCESSING  A  PICTURE  IF  AN  ERROR  IS  DETECTED 

•SUCH  AS  WHEN  A  ROW  OR  COLUMN  IS  TO  BE  PROCESSED  -  LESS  THAN  0  OR  GREATER  THAN 
256.  THE  LIMITS  ARE  AUT OMAT I C AL L V  CHANGED  TO  THE  MlN/MAX  ALLOWABLE  THIS 
PROCEDURE  ALSO  ADJUSTS  THE  LIMITS  BASED  ON  THE  SIZE  OF  MASK  USED  AND 
MIN  MAX  EDGE  BOUNDARIES  PERMITTED  F OR  PROCESSING  *) 

PROCEDURE  CHE'CKLIMITS  (VAR  ORCW.CKCOL  ROWCOL); 


BEGIN 


IF 

MASK  = 

0 

THEN  SIZE 

=  0  ELSE 

IF 

MASK  = 

3 

THEN  SIZE 

-  1  ELSE 

IF 

MASK  = 

5 

THFN  SIZE 

=  2  ELSE 

IP 

MASK  =* 

7 

THEN  SIZE 

-  3  ELSE 

IF 

MASK  =* 

9 

THEN  SIZE  : 

=  A  ELSE 

IF 

MASK  » 

1  1 

THEN  SI ZE 

=  5. 

IF  (  ( C  k  ROW 
PE  C  IN 

WHILE 

~ 

SIZE'  C  0)  THEN 

( < CKRQW  -SIZE  )  <  < 

BEGIN 

CKROW  =  C*ROW  *  1, 

END. 

Wfi;TELN<-  CKROW  LIMIT  7  JO  LOW  -  CHANGED  TO  '.CKROW), 
M INROW  =  CKPOW. 

END. 


IF  MC*»DW  ♦  SIZE)  J  255)  T^EN 
BEGIN 

•  'C.MP2W  -  SIZE  )  D  * 
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/ 


3EGIN 

CKRQU  *  CKROW  -  1. 

END. 

WR I TELN ( '  CKROW  LIMIT  TOO  HIGH  -  CHANGED  TO  CKROW), 
MAXROW  =  CKROW. 

END. 


IF  ((CKCOL  -  SIZE  )  <  0)  THEN 
BEGIN 

WHILE  <  (  CKCOL  -  SIZE)  <  0  )  DO 
SEGIN 

CKCOL  =  CKCOl  ♦  1. 

END. 

WRJTELNt'  CKCOL  LIMIT  TOO  LOW  -  CHANGED  TO  '.CKCOL)* 
MINCOL  =  CKCOL. 

END. 

IF  ((CKCOL  ♦  SIZE  >  >  255  )  THEN 
BEGIN 

WHILE  ((CKCOL  ♦  SIZE)  >  255)  DO 
BEGIN 

CKCOL  ■*  CKCOL  -  l. 

END* 

WR ITELN( 'CKCOL  LIMIT  TOO  HIGH  -  CHANGED  TO  '.CKCOL); 
MAX  CO’.  *  CKCOL. 

EJ '  j 


IF  ANSWER  = 
BEGIN 
M INROW 
MAXROW 
MINCOL 
MAXCOL 
END. 


YESCHAR  THEN 


-  SIZE. 

=  255  -  SIZE. 
=  SIZE. 

=  255  ~  SIZE. 


END. 


PROCEDURE  SETUP. 
LABEL 

10.  20  . 


VAR 

DUM1.DUM2  INTEGER. 


BEGIN 

DUM 1  =  1. 

DUM2  *  I; 

WRIT  ELN. 

WRJTELNt  'DO  YOU  WANT  TO  FILL  THE  OUTPUT  DISK  FILES  -'), 

WR  I  TELN  (  '  NEWD1  TO  NEWD4  WITH  BLANKS’7  Y/N  '); 

READLN( ANSWER), 

IF  ANSWER  =  \ ESCHAR  THEN 
BEGIN 

WR  J  TELN, 

WR  I  1  ELN  (  '  THE  "WORKi^G,‘  DISK  FILES  ARE  NOW  BEING  FILLED  WITH  BLANKS 
BLANKFILE. 

WR  I  "ELN. 

WRITELNt '  NEWD1.NEWD2  -  READY  '), 

BLANK2F 1LE. 

WR I ” ELN ( '  NEWD3.NEWD4  -  READY'). 
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w»I  7ElN. 


END. 

10  WRITELNt'  ALLOWABLE  mask  S17ES  F OR  PICTURE  PROCESSING  ARE  3X3.  5X5.  7x7.  9x9. 

WRITE*  ENTER  MASK  SIZE  ,,3“  “V  "7"  "9"  "  1  1  “  '). 

READLN  i  MASK) . 

IF  <  MASK  O  3)  AND  (MASK  O  5)  AND  (MASK  O  7)  AND  (MASK  <>  9)  AND  (MASK  <>  Jl> 
BEGIN 

WR  ITELN, 

GOTO  10. 

END. 

WnlTELNC  MASK  SIZE  ’.MASK.  '  BFING  USED'). 

-RITfLN. 

20  WRIT£LN(  '  DO  YOU  WANT  TO  PROCESS  ENTIRE  PICTURE  -  Y  OR  N  *■>'). 

RE  ADEN ( ANSWER > . 

IF  ANSWER  =  NOCHAR  THEN 
BEGIN 

WR ITELN. 

WRITFLNt'  ENTER  -  "  MINROW  MAXROW  MINCOL  MAXCQL  '); 

'WRITELNt'  MIN/ HA X  LIMITS  ON  VALUES  ENTERED  *  0  THRU  255  INCLUSIVE'). 

WRITELNt'  FOR  EXAMPLE  2  20  3  125  '). 

R E A DL N ( M I  NR OW ,  MAX ROW.  M I NC OL .  HA X C OL ) . 

CMECKLIHITS<M2NROW.  MINCOL). 

WR ITELN. 

CHECKLIMI7S< MAXROW,  MAXCOL). 

WRITELNt  MINROW  MINROW,  '  MINCOL  '.MINCOL). 

URIT£LN(  MAXRQW  '.MAXROW.  '  MAXCOL  '.MAXCOL). 

WR  I  TELN, 

END 

ELSE 

IF  ANSWER  =  YfSCHAR  THEN 
DEGIN 

CHECKL I M I T  5 ( DUM 1 ,  DUM2 ) . 

IF  MASK  O  3  THEN 

WRITELNt'  IGNORE  ERROR  MESSAGES  ON  ROW/COL  P ARAMET  ERS ' ) ; 

WR I TELN< '  PICTURE  WILL  PROCESS  FROM  ROW  '.MINROW,  'TO  '.MAXROW); 

WR I T  ELN  <  '  AND  FROM  COL  ".MlwCOL.  '  TO  ' .  MAXCOL ) ; 

END 
Ft  SE 
DEC  I N 

WR I TELN.  WR ITELN.  WR ITELN.  WR I TELN; 

WR  I  TEl.N  (  '  INCORRECT  ENTRY  -  REPLY  WITH  "V  OR  “N"  -  TRY  AGAIN1'). 

GOTO  20. 

END. 

WRI TELN. 

WRITELNt*  SETUP  COMPLETE  «****'>, 

WRI TELN. 

END. 


<*  PROCEDURE  THRESHOLD  IS  DESIGNED  TO  FIND  BOTH  THE  MEAN  AND  57ANDARD  DEVIATION 
OF  ALL  THE  PIXELS  IN  A  256  BY  256  PICTURE  AT  THIS  TIME  ALLOWANCE  HAS 
NUT  BfEN  MADE  FOR  DOING  SUB  IMAGE  AVERAGING.  SUCH  AS  MAY  BE  SPECIFIED  BY 
THE  LIMITS  MIN.  MAXROW  AND  MIN.  MAXCOL  Bt CAUSE  OF  INACCURACIES  ON  THE  LEFT 
AND  BOTTOM  BORDERS  Th£  AVEPACE  IS  COMPUTED  FROM  12  TO  256  AND  FROM 
O  TO  233  *) 


OP  11X11 

THEN 


i. u, r. ;  r ThKc  5H0LD# 


VAR 

THRFILE  DATAFILE. 

TS3R.  T SUM.  TNUM.  P I XPT  REAL. 

T1F1LE. T2FILE. T3FILE. TAFILE  DAT AF . 

TAVC. STANDEV  RFALi 

T BLOCK  PACKED  ARRAY  TO  633  OF  DATAFILE. 

THRROU.  T  HR  COL  O  255. 


'  TSQR  =  O.  <»  EQUALS  THE  sun  OF  the  pixel  values  squared* I 

PIXPT  =  0.  (.COUNTS  the  »  OF  PIXELS  PROCESSED*) 

isun  *  0.  (.EQUALS  sun  OF  ALL  PIXEL  VALUES*) 

F IL£TITLE< T1FILE. 'DATA1'). 

RESET  <  7  IF  1 LE ) » 

FOR  THRROW  5  0  TO  63  DO 
BEGIN 

get<tifile>. 

TBLOCKC THRR0W3  *  T1FILE''. 

THRFILE  =  TBLQCKC  THRROW 3  i 
WITH  THRFILE  DO 


BEGIN 


FOR  THRCOL 
BEGIN 

TNUM 

tsum 

P  I  XPT 

T5QR 


=  13  TO  255  DO 

*  LINEC  THRCOL  3 i 
»  TSUM  ♦  TNUM; 

a  PIXPT  ♦  1, 

=  T5GR  ♦  < TNUM  • 


TNUM). 


END. 

END. 

CLOSE(TIFILE) . 

WR  I  TELN. 

UR 1TELN<  '  1 '4  DONE  ’  ) . 

FILETITtEt T2FILF.  ’  DATA2  '  ) . 

RESET ( T2FILE) - 

FOR  THRROW  =  0  TO  63  DO 

ELGIN 

GET i 7  ZF I LE  > ; 

TBLOCKC THRROWT  =  T2FILE- 
THRFILE  =  TBLOCKt THRROWT. 

WITH  THRFILE  DO 
BEGIN 

FOR  THRCOL  ■ *  13  TO  255  DO 
BEGIN 

TNUM  =  LI  NEC THRCOL 3. 

TPUM  =  TSUM  ♦  TNUM. 

PIXPT  =  PIXPT  1. 

TSGR  =  TSQR  ♦  <  TNUM  *  TNUM). 

END 

END. 


END. 

Close < tzfile). 

wR I TELN ( '  1/2  DONE  >; 

F IlET I TLE ( T3FILE.  'DATA3 ' ) » 

RL5E  T  <  T3F ILE ) * 

rfiR  '"HRROW  *  O  TO  63  DO 

BEGIN 

GET/T3FILE). 

TBLOCKt  THRROW]  *  T  3F  ILE 
THRFILE  *  TBLOCKl THRR0W3. 
WITH  THRFILE  DO 
BEGIN 

FOR  THRCOL  -  13  TO  255 


DO 
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BEGIN 

TNUM  =  L I  NEC ThrCOLJ. 

IS UH  *  T  SUM  *  TNUM. 
p I XPT  =  P  I  XPT  ♦  1 . 

T5QR  *  T50R  ♦  C  TNUM  *  TNUM  ). 

END 

END. 

END. 

Close  <  T3F i  le  ) . 

WR 1TELN( '  3/4  DONE'). 

UR  ITELN. 

FILE TITLE ( T4FILE.  ' DAT A4 ' ) . 

Pi  SE  T  <  T4F I LE  )  i 

FOR  TwRROW  *  0  TD  40  DO 

BEGIN 

GET  <  T  4F I LE ) . 

TBLOCKt  7 HR ROW 3  =  T4F I LE  ^ . 

THRFILE  *  TBLOCKt TnRROW). 

WITH  THRFILE  DO 
BEGIN 

FOR  THRCOL  *  J  3  TO  D55  DO 
BEGIN 

TNUM  »  L INEI ThRCOL 3. 

TSUM  «  TSUM  ♦  TNUM. 

P I XPT  =  P  I  XPT  ♦  1 . 

TSQR  «  TSQR  ♦  (  TNUM  *  T NUM ) . 

END 

END. 

END. 

CLOSE \T4FJLEJ» 

TAVG  *  <  T SUM  /  PI  XPT). 

MEAN  =  ROUND (TAVG). 

WR  I  TELN<  '  MEAN  PIXEL  VALUE  F- OR  256  X  25e>  PICTURE  =  '  .  T  SUM.  '  /  '  .  P  I  XPT. 

< *CuMPUTE  standard  DEVIATION*) 


STANDEV  =  SQRT<  < (PiXpT*TSQR)  -  < 7  SUM* 7  SUM ) ) / ( P I X P T * < P J XP T - \  )  >  ) , 
STDEVIATION  =  ROUND< STANDEV)  , 

WRI 7ELN<  '  57 ANPARD  DEVIATION  =  STDEVIATION). 

END. 


MEAN)  , 


) 


<*  PROCEDURE  READPIX  IS  DESIGNED  Tq  RETURN  A  GRFVlFVFL  CAiut  F DR 
EACH  PIXEL  IN  A  PICTURE  MATRIX  OF  256  D y  256  DISCRETE  PICTURE 
p0 1  NTS  THE  PARAMETERS  PASSED  INTO  ThF  PROCEDURE  ARE  THE  X.Y' 
COORDINATES  OF  THE  POINT/PIXEL  IN  THE  DIGITIZED  PICTURE*) 

PROCEDURE  RE ADP I  X (  P1XR0W  ROUCOL, 

VAR  ROWDATA  RECROWDAT  A ) i 

(*  °0UDAT A  IS  THE  DATA  WHICH  IS  PASSED  SACK  FROM  READPIX  IT 
CONSISTS  OF  A  PACKED  RECORD  OF  PICTURE  ROWS  -  EACH  CONTAINING  256 
DISCRETE  PICTURE  GPEYLEVEL  VALUES  THE  SIZE  QF  ROWDATA  IS  COMPLETELY 
DEPENDENT  ON  THE  "MASK "  SIZE  EITHER  3X3  -  3R0WS.  5X5  -  5  ROWS. 

OR  ?X7  -  7  ROUS  *) 


.  AR 

DONE,  DTUO.  DTHREE,  OF OUR  DA T AF . 

RECINDEX  INTEGER. 

NUMPI XffOW.  70PR0W. 80TRDW  INTEGER. 


SEC  INDEX  =  0. 

IF  PIXROW  <  0  THEN 

wRl'ELN  <  'ERROR  JN  PROCEDURE  READPJX  ~  VALUE  <  1'); 
IF  PIXROW  D  255  THEN 

WRJTFLNt'  ERROR  IN  PROCEDURE  READPIX  -  VALUE  >  255') 


IF  MASK  »  0  THEN  OFFSET  »  0  ELSE 

IF  MASK  =  3  THEN  OFFSET  •  =  1  ELSE 

IF  MASK  »  5  THEN  OFFSET  *  2  ELSE 

IF  MASK  »  7  THEN  OFFSET  -  3  ELSE 

IF  MASK  *  9  THEN  OFFSET  :  *  4  ElSE 

IF  MASK  *  11  THEN  OFFSET  *  5. 

TOPROW  =*(  PIXROW  -  OFFSET). 

BOTROW  »<  PIXROW  ♦  OFFSET ) ; 

SIZE  *  OFFSET, 

FDR  PIXROW  »  TOPROW  TO  BOTROW  DO 
BEGIN 

REC INDEX  =  REC1NDEX  ♦  1; 

IF  (PIXROW  AND  (PIXROW  <=  63)  THEN 

BEGIN 

FILET ITLE (DONE.  'DATA1 ' ) ; 

RESET (DONE), 

SEEK (DONE, PIXROW); 

GET (DONE); 

ROWDATAC  REC  1 NDEX  3  .  *  DONE''. 

CLOSE ( DONE ) ; 

END 

ELSE 

IF  (PIXROW  >=*  64)  AND  (PIXROW  127)  THEN 
BEGIN 

NUMPIXROW  =  PIXROW  -  64; 
FILETITLE(DTW0.  '  DATA2  ' ) . 

RESET ( DTWO ) i 

SEEK ( DTWO, NUMPIXROW), 

GET (DTWO). 

ROWDATACREC INDEXI  »  DTWO~; 

CLOSE ( DTWO ) ; 

END 

ELSE 

I F ( P 1 XROW  >*  128)  AND  (PIXROW  <*  191)  THEN 
BEGIN 

NUMPIXROW  «  PIXROW  -  12B. 
FILETITLE(DTHREE.  ' DAT  A3  '  ) ; 

RESET ( DTHREE ) ; 

SEEK ( DTHREE, NUMPIXROW); 

GET (DTHREE); 

RQWDAT A [ REC I NDEX )  »  DTHREE^. 

CLOSE ( DTHREE ) ; 

END 


IF  (PIXROW  19?)  AND  (PIXROW  <»  255)  THEN 
BEGIN 

NUMPIXROW  «  PIXROW  -  192. 

F I  LET  I TLE ( DFOUR .  DAT  A4  '  ) . 

RESET ( DF OUR ) , 

Sc  Eh  I DFOUR.  NUMP I XROW) . 

GE  T ( DF  OUR  ) . 

ROUCATAt REC INDEX!  »  L  OUF'. 

Ct-LSf  ( l/FOUB  >  • 
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END. 

END. 

END. (♦READPlX#) 

<  •»««««  ********  mmm  »  mm  ******  «*««*»•«#  »»  »  »«»»  «»  *«  *«  •«  «  ««  * «  ««  ««  ««  «*«  ««««  «  #«####•#».#•«#♦  ) 

(♦PROCEDURE  STORELINE  IS  DESIGNED  TO  STORE  ONE  COMPLETE  LINE  (OR  ROW)  OF 
A  PICTURE  (IE  256  PIXELS)  INTO  A  DISK  FILE  THE 
DISK  FILE  IS  COMPOSED  OF  64  SUCH  LINES  AND  IS  OF  TYPE  DAT AF  AS 
DESCRIBED  ABOVE  EACH  LINE  .ITSELF,  IS  A  PACKED  RECORD  OF  PIXELS  OF 
TYFE  ■,DATAFILEm  .  *> 


PROCEDURE  STORELINE (ROWSTORE  ROWCOL. 

DATALINE  DATAFILE). 


VAR 

DIFILE.  D2FILE. D3FILE, D4FILE  DATAF. 
STROBJECT  :  DATAFILE; 


BEGIN 

I F ( ROWSTORE  >*0)  AND  (ROWSTORE  <=  S3)  THEN 
BEGIN 

STROBJECT  . *  DATALINE. 

FILETITLE(D1FILE.  'NEUD1  ' ); 

RESET (DIFILE); 

SEEK i DJFILE.  ROWSTORE); 

WRITE  (DIFILE.  S'lROBJECT); 

CLOSE  (DIFILE); 

END 

ELSE 

IF ( ROWSTORE  >=  64)  AND  (ROWSTORE  <=  127)  THEN 
BEGIN 

ROWSTORE  =  ROWSTORE  -  64; 

STROBJECT  .=  DATALINE; 

FrLETITLE(D2FlLE.  'NEWD2 ' ) . 

RESET ( D2F I LE ) ; 

SEEK(D2FILE. ROWSTORE). 

WR  2  TE ( D2F ILE. STROBJECT). 

CLOSE (D2FILE). 

END 

ELSE 

1 F ( ROWSTORE  >«  120)  AND  (ROWSTORE  <=  191)  THEN 
BEGIN 

ROWSTORE  =  ROWSTORE  -  128, 

STROBJECT  ;«  DATALINE, 

FILETITLE(D3FILE.  NEWD3'); 

RESET  <  D3FILE  > ; 

SEEK ( D3FILE, ROWSTORE); 

WR I TE ( D3FILE. STROBJECT); 

CLOSE  (D3FILE); 

END 

ELSE 

IF (ROWSTORE  >*  192)  AND  (ROWSTORE  <-  255)  THEN 
BEGIN 

ROWSTORE  : *  ROWSTORE  -  192, 

STROBJECT  =  DATALINE. 

FILET1TL.E(D4FILE.  NEWD4'). 

RESET(D4FILE) . 

SEEK  <  D4F ILE. ROWSTORE). 

WR  1TE(D4F  ILE.  STROOv)ECT  ) , 

CLOSE  <  D4F I LE ) • 

END 


i 


CL  fit 

*RITELN(  "'ROWSTDRE''  VALUE  PASSED - OUT  Ob  RANGE  IN  ST  OR  EL  1 NE  ’  )  • 

END- 


1  . nn*mnn*m** ******* **H ********  ****************** •***•*•*•****•»* »**•*»»*••*******  > 


'♦PROCEDURE  FILLMAS*  TAKES  A  SPECIFIC  PIXEL  IN  THE  PICTURE  AND  BUILDS 
a  MATRIX  OR  'MASK'  WITH  THE  SPECIFIED  PI  'L  AS  THE  CENTER 
MATRIX. MASK  SIZE  IS  ENTERED  INTERACTIVELY  BY  THE  USER 
ALLOWABLE  MASK  SIZES  ARE  3X3,  5X5,  7X7.  9X9  OR  11X11 

THIS  PROCEDURE  DOES  NOT  COMPENSATE  FOR  ’EDGE'  PIXEL  VALUES* ) 

PROCEDURE  FILLMASK<  FROW  ROWCOL. 

PIXPOINTS  :  RECROWDATA. 

VAR  FMASK  :  RECMATR I XMASK ) i 

V  AR 

FCOL,  ECOL, SCOL  .  ROWCOL; 

TOPCORNER,  botcorner. leftcol,  RTCOL,  LEFTCORNER  INTEGER. 

F INDEX,  FCOUNTER.  FNUMBER  :  INTEGER; 

FILLPIXREC  0  .  255; 

FSIZE. MARKER  INTEGER, 

TEST INDEX  1  1 21; 

BEGIN 

FSIZE  :«  OFFSET; 

FCOL  : *  MINCOL, 

SCOL  =  MINCOL. 

ECOL  =  MAXC0L; 

TOPCORNER  : =  FROW  -  FSIZE. 

BOTCORNER  •=  FROW  +  FSIZE; 

LEFTCOL  =  MASK  -  FSIZE; 

RTCOL  *  MASK  *  FSIZE; 

LEFTCORNER  =  FCOL  -  SIZE. 

MARKER  : -  LEFTCORNER, 

<*  THIS  PROCEDURE-  LABELS  the  POINTS  in  THL  MATRIX  /  MASK 
STARTING  AT  ONE  "l"  UP  TO  THE  (MASK  ♦  MASK)  LOCATION 
NOTE  HOWEVER  THAT  RECORDS  OF  THESE  MATRIX  FILES  GO  FROM 
O  TO  255  *) 

BEGIN 

FOR  FILLPIXREC  =  SCOL  TO  FCOL  DO 
BEGIN 

F1NDEX  -  O; 

FCOUNTER  »  O; 

FOR  FROW  . =  TOPCORNER  TO  BOTCORNER  DO 
BEGIN 

FCOUNTER  : *  FCOUNTER  ♦  1; 

F NUMBER  : -  MARKER; 

FOR  FCOL  -  LEFTCOL  TO  RTCOL  DO 
BEGIN 

FINDEX  : *  FINDEX  ♦  1; 

FMASK CF ILL P I XREC  3  P I XMATR I  X  IF  INDEX  3  «  P I XPO I  NTSC FCOUNTER 3  LINE! FNUMBER 3 • 

F NUMBER  FNUMBER  *  1. 

END 

END; 

MARKER  =  MARKER  ♦  1; 

END. 

END; 

<*  FOR  FILLPIXREC  *  MINCOL  TO  maXCOL  DO 

BEGIN 

WRITELN< 'FMASK  USING  MASK  '.FILLPIXREC.  '  AND  MASK  SI 7E  ’.MASK). 
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1  TO  (MASK  *  MASK)  DO 


FGR  TEST  INDEX  = 

BEGIN 

wR I T£ * FMASK  CF1LLP1XREC]  P I XMATR I  XI  TEST  INDEX  ]  3), 

END. 

URITELN. 

URITELN*  '.##«#*♦*#**♦*♦♦*♦♦*•#**♦♦****♦*♦**♦*.••*♦*♦♦#♦*' ) , 

END.  *) 

END. 


i  •««««*»«  «*«*  mm*  «*  »««  •«««  ######*  *>  a  #  *  ««»  mm  mm  *  mmm  m  m  mm  *  m  *  mm  m  m  mm  m  m  mm  m  mm  m  mmm  mmmi-  mm  m  *  m  m  *  *  m «  mm  m  m  m  mm  m  m  m  m  m  *  > 


PROCEDURE  AVERAGE (AVGMASK  REC MATR 1 XMASK ) . 


AVGFILE  DATAFILE. 

AVGREC  0  255  ; 

AVERAGE.  AVDEX  INTEGER. 

AVCSUM  REAL. 

AVP  PACKED  ARRAYC I.  .  1213  OF  0 
MASKSQUARED  :  INTEGER. 


BEGIN 


WRITELN*  PROCESSING  WITH  A  M.  SK  SIZE  OF  '.MASK); 

AVERAGE  =0. 

MASKSQUARED  :  =  MASK  *  MASK, 

FOR  AVGREC  *  MINCOL  TO  MAXCOL  DO 
BEGIN 

AVCSUM  =  O. 

FOR  AVDEX  1  1  TO  MASKSQUARED  DO 
BEGIN 

AVPCAVDEXJ  =  AVGMASK C AVGREC I  P I XMATR I  X C AVDEX 3 . 
AVCSUM  *  AVGSUM  +  AVPt AVDEX]; 

END. 

AVERAGE  =  ROUND ( AVGSUM/ MASKSQUARED )  ; 

.  GF1LE  LINECAVGREC3  *  AVERAGE, 

END, 

STCRELINE^ROW.  AVGFILE), 

AVEFAGE  =  0; 


END, 

<  *## 


BEGIN  <*  MAINLINE*) 

LINECOUNT  :■  Oi 

WRITELN.  URITELN.  WRITELN, WRITELN.  WRIT  ELN;  WRITELN.  URITELN,  WRITELN, 

WR I TELN<  '  WELCOME  TO  PROCESS  PICTURE  -  YOU  ARE  CURRENTLY  USING  THE  PASCAL  PROGRAM.  .  . 
WRITELN; 

UR  I TELN  *  '  " AVGMASKS"  (VERSION  MAY  1982)'); 

URITELN. 

WRITELN. 

WRI TELN (  'NOTE  FOR  THIS  PROGRAM  TO  WORK  A  PICTURE  FOR  PROCESSING') 

WRITELN<  'MUST  ALREADY  EXIST  ON  DISK  IN  FILES  SPECIFICALY  LABELLED  -  DATA1, 

WR I TELN (  '  DATA?,  DAT  A3  AND  DATA4 ' )  . 

WRJTFLNf'TO  CREATE  THESE  FILES  SIMPLY  ENVOKE  THE  FORTRAN  PROGRAM  "V4"  BY  TYPING  "V4" 
WR I TELN ( 'AT  THE  TERMINAL.  YOU  MUST  EXIT  THIS  PROGRAM  FIRST,  HOWEVER.  SO ' ) I 
WR I TELN  <  '  TYPE  "CONTROL  A",  AND  REENTER  THIS  PROGRAM  AFTER  DATA I  THRU  4  ARE  CREATED' 
URITELN.  WRITELN.  WR I TELN,  WRITELN,  WRITELN, 

82  WRITELN* 'BESIDES  PERFORMING  LOW  PASS  (AVERAGING)  FILTERING  OPERATIONS  '). 

WRITELN* 'THIS  PROGRAM  CAN  CALCULATE  THE  MEAN  AND  STANDARD  DEVIATION'); 

URITELN*  'FOR  THE  PIXEL  VALUES  OF  THE  ENTIRE  PICTURE.  DO  YOU  WISH  TO'); 

WRITE*  'HAVE  THESE  VALUES  CALCULATED?  Y/N  '). 
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readln-  answer ) . 

IF  ANSWER  =  YESCmAR  Then 
BEGIN 

WRITELN. 

WRITELN!  '  ROGER.  USER  WAIT  ONE  '). 

THRE5HOLD, 

WRITELN.  WRITELN.  WRITELN.  WRITELN,  WRITELN.  . 

WRJTELNf'  PHEW  *#«•  .  -  I  HOPE  I  DO  NOT  HAVE  TO  DO  THAT  AGAIN  ' 

WRJTELN. WRITELN. WRITELN, WRITELN. WRITELN, 

END 

else 

IF  ANSWER  O  NOCHAR  THEN  GOTO  82. 

WRITELN; 

SETUP, 

FOR  ROW  •*  MINROW  TO  MAXROW  DO 
BEGIN 

WRITE! 'ENTERING  READPIX  -  ROW  *  '.ROW); 

READP IX! ROW,  MATRIXDATA); 

WRITELN!'  ENTERING  FILLMASK  -  ROW  »  '.ROW); 

FILLMASK ( ROW,  MATR I XDATA.  GRECMASK ) ; 

AVERAGE ! GRECMASK ) ; 

WRITELN!  'LINECOUNT  -  ' .  L I NECOUNT > . 

LINECOUNT  : *  LINECOUNT  >J. 

END; 

WRITELN.  WRITELN; 

WRITELN!'  PROCESSING  OF  PROGRAM  AVCMASKS  COMPLETE 

WRITELN. WRITELN; 

WRITELN!'  THANK  YOU.  HAVE  A  GOOD  DAY*'). 

WRITELN. WRITELN, 

END 
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PROGRAM  name  PROCESS  PICTURE  "AN DOR  PA" 

CAPT  BLAINE  FELTMATE 

DATE  SEP  1981  REVISED  21  MAY  1902 

PURPOSE  THE  PURPOSE  OF  THIS  PROGRAM.  ANDQR  PA.  IS  TO  PERFORM  BOOLEAN 

"AND"  AND  "OR”  OPERATIONS  ON  TWO  INPUT  IMAGE  FILES  AND  CREATE  A  THIRD 
OUTPUT  IMAGE  FILE  OF  THE  RESULTING  PROCESSED  PICTURES  THE  FIRST  INPUT  IMAGE 
FILE  TO  THE  ANDOR  PROGRAM,  SIMPLY  CALLED  IMAGE  ONE. 

IS  CREATED  BY  RUNNING  THE  FORTRAN  PROGRAM  VP  I C  BY  TYPING  "VPIC"  AT  THE 

terminal  this  must  be  done  before  entering  this  program  vpic  takes  a 

64K  PICTURE  FILE  AND  DIVIDES  IT  EVENLY  INTO  FOUR  16K  FILES  WHICH 
CAN  BE  USED  FOR  PROCESSING  BY  THIS  PROGRAM  THIS  IS  DUE  TO  THE  l&K  STACK 
limitation  IMPOSED  BY  THE  RDOS  VERSION 

of  pascal  running  on  the  afit  signal  processing  laboratory  computer 
vpic  then  creates  files  pici  thru  pica  qn  disk  the  SECOND  INPUT 

IMAGE.  IMAGE  TWO.  IS  PREPARED  FOR  PROCESSING  IN  A  SIMILAR  MANNER. 

TC  CREATE  the  a  disk  files  for  IMAGE  TWO  SIMPLY  TYPE  ,,VD"  and  files 
Dl  THRU  D4  WILL  BE  CPEATED  AGAIN  THIS  5HOULD  HAVE  BEEN  DONE  PRIOR  TO 
ENTERING  THIS  PROGRAM 

THE  OUTPUT  OF  THE  ANDOR  PROGRAM.  EITHER  "AND"  OR  "DR" 

OPERATION  SPECIFIED  BY  THE  USER.  WILL  BE  STORED  IN  FILES  Nl.  N2.  N3  AND 
N4  ON  DISK.  SIMPLY  TYPE  "APPEND  <FILENAME>  Nl  N2  N3  N4  ■' 

AT  THE  TERMINAL  TO  REFORM  THE  4  PROCESSING  FILES  BACK  INTO  A  SINGLE 
PICTURE  FILE  ALL  OTHER  PROMPTS  SHOULD  BE  SELF  EXPLANATORY  AS  THEY 
APPEAR  INTERACTIVELY  TO  THE  USER 

NOTE  THE  VALUES  STORED  IN  THE  OUTPUT  FILES  ARE  TAKEN  FROM 
THE  INPUT  IMAGE  TWO  FILES  < VD  -  Dl  D2  D3  D41 


*«*«««*««  *«««««*#**»#  **««*»«  *#*»**«**«*  «*«**»*«««#*«  ««#<**««**«  *«***##«#»#»*»**»»*»> 


PROGRAM  P  I  X VALUE  < INPUT.  OUTPUT.  DATA1.  DATA2.  DATA3. DA7A4, 
Dl. D2- D3, D4. P IC1. PIC2. PIC3. PICA. 

Nl . N2, N3. N4 ) i 


LABEL 

1111; 

CONST 

VESCMAR  »  'Y ' , 

NOCHAR  *  'N 
ANDCHAR  *  'A'# 

ORCHAR  *  0  '# 

type 

DATAFILE  *  PACKED  RECORD 

LINE  .  PACKED  ARRAYCO.  .  2553  OF  0.  .  IS 
END; 

DATAF  «  FILE  OF  DATAFILE; 

ROWCOL  *  0  255; 

^AR 

DATAl .  DATA?.  DATa3,  DATA4  DATAF, 

Nl,  N2.  N3.  N4.  Dl.  D2.  D3.D4.P1C1.P1C2.PIC3.PIC4  DATAF; 
ANSWER  CHAR; 

MEAN, STDSVIATION  INTEGER; 

ANDTYPE  BOOLEAN; 

KONST  REAL; 
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<*  PROCEDURE  THRESHOLD  JS  DESIGNED  TO  FIND  BOTH  THE  MEAN  AND  STANDARD  DEVIATION 
OF  ALL  THE  PIXELS  IN  A  256  BY  256  PICTURE  AT  THIS  TIME  ALLOWANCE  HAS 
NOT  BEEN  MADE  FOR  DOING  SUB- IMAGE  AVERAGING.  SUCH  AS  MAY  BE  SPECIFIED  BY 
THE  LIMITS  MIN,  MAXROW  AND  MJN,  MAXCOL  BECAUSE  OF  INACCURACIES  ON  THE  LEFT 
AND  BOTTOM  BORDERS  THE  AVERAGE  IS  COMPUTED  FROM  12  TO  256  AND  FROM 
0  TO  233  *> 


procedure  threshold. 


VAR 

thrfile  datafile. 

PJXPT,  TNUM,  TSUM,  TSGR  REAL.  (^REQUIRED  DUE  TO  SIZE  OF  NUMBERS  WHICH  ACCUMULATE* 1 
TlFlLE. T2FILE, T3FILE,  T4FILE  DATAF. 

TAVG. STANDEV  REAL, 

tblock  .  packed  array  ro  633  of  datafile* 

THRROW. THRCOL  .  0  255* 

BEGIN. 

TSQR  =*  Oi  (*  EQUALS  THE  SUM  OF  THE  PIXEL  VALUES  SQUARED* ) 

PIXPT  *  O*  (#COUNTS  THE  «  OF  PIXELS  PROCESSED*) 

TSUM  *  0.  <  *E DUALS  SUM  OF  ALL  PIXEL  VALUES*) 

FILETJTLE<71FJLE.  ' DATA1  ' ) » 

RESET ( T IF  ILE ) * 

FDR  THRROW  *  0  TO  63  DO 
3EGIN 

GET(TIFILE). 

TBLOCK  I  THRROW!  =  TJFILE*'* 

THRFILE  *  TBLOCKC THRROW!, 

WITH  THRFILE  DO 
BEGIN 

FDR  THRCOL  =  13  TO  255  DO 
BEGIN 

TNUM  =  LI NEC THRCOL 3. 

TSUM  =  TSUM  *  TNUM, 

PIXPT  *  PIXPT  1. 

TSQR  .  *  TSQR  ♦  (TNUM  *  TNUM). 

END, 

END, 

END* 

CLOSE  <  TJF JLE ) , 

WR I TELN< '  1/4  DONE'). 

FILETITLEC T2FILE.  DATA2'); 

RESET \T2F1LE)( 

FOR  THRROW  .  *  O  TO  63  DO 
BEGIN 

GET ( T2F  TLE ) . 

TBLOCKC  THRROW 3  »  T2FJLE', 

THRFILE  3  TBLOCKC THRR0W3* 

WITH  THRFILE  DO 
BEGIN 

FOR  THRCOL  «  13  TO  255  DO 
BEGIN 

TNUM  -  LINE! THRCOL 3, 

TSUM  *  TSUM  ♦  TNUM. 

PIXPT  «  PIXPT  ♦  1, 

TSQR  *  TSQR  ♦  (  TNUM  «  TNUM). 

END 

END, 

END. 
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.  .  -  -  ,  "V'  * 


C;.CSc  »  72F JLE)  . 

wR!TELN(  '  1/2  DONE  '  ). 

FlLETJTLEt T3FILE.  '  DATA3 ' ) » 

RESET ( T3F I LE ) ; 

FOR  T HR ROW  =  O  TO  63  DO 
BEGIN 


GET ( T3F I LE ) . 

TBLOCKCTHRROW3  =  T3FILE~; 

THRFILE  =  TBLOCK [ THRROW]  , 

WITH  THRFILE  DO 
BEGIN 

FOR  THRCOL  :  -  13  TO  255  DO 
BEGIN 

TNUM  :■  LI NEC THRCOL 3; 

TSUM  :*  TSUM  ♦  TNUM. 

PIXPT  : =  PIXPT  ♦  li 

TSQR  :  =  T SGR  ♦  (  TNUM  *  TNUM  >; 

END 

END. 

END; 

CLOSE  <  T3F I LE ) ; 

WR I TELN< '  3/4  DONE'). 

FILETITLE(T4FILE.  ' DATA4 ' ) , 

RESET  <  T4FILE  ) . 

FOR  THRROW  =  O  TO  40  DO 
EECIN 


GET ( T 4F ILE  ) . 

TBLOCK  C  THRROW 3  :  =  T  4F  2 LE 
THRFILE  ®  TBLOCK  C  THRROW 3 ; 

WITH  THRFILE  DO 
BEGIN 

FOR  THRCOL  =  13  TO  255  DO 
BEGIN 

TNUM  *  LINED THRCOL 3; 

TSUM  .  *  TSUM  ♦  TNUM, 

PIXPT  =  PIXPT  -f  1; 

T  SGR  *  T  5QR  ♦  (  TNUM  *  TNUM). 

END 


FND. 

END. 

CLC^E  t  T4F I L  E ) , 

WRITELN( '  4/4  DONE ' ); 

TAVC  =  (TSUM  /  PIXPT). 

MEAN  =  ROUND (TAVC). 

WR I TELN <  '  MEAN  PlXELVALUE  FOR  256  X  256  PICTURE 


(♦COMPUTE  STANDARD  DEVIATION#) 


'.  TSUM. 


/  PIXPT.  ' 


5TANDEV  =  SORT ( ( (PIXPT#TSQR)  -  ( TSUM#T 5UM ) ) / < P I XPT ♦ ( P I XPT - 1 ) ) ) . 
5TDEVI AT  I  ON  *  ROUND ( ST ANDEV ) ; 

UR  I TELN (  '  STANDARD  DEVIATION  =  ' .  S TDEV I  AT  I  ON ) ; 

WR I TELNj 

END; 


(♦PROCEDURE  MATCHPIC  DOES  THE  ACTUAL  BOOLEAN  ••AND”  AND  "OR"  OPERATIONS 
THE  CHOICE  OF  WHICH  OPERATION  IS  PERFORMED  IS  ENTERED  AT  THE  TERMINAL 
B V  THE  USER.  AFTER  AN  APPROPRIATE  PROMPT  REQUESTING  THIS  INFORMATION 
IS  GIVEN  BY  THE  COMPUTER  ♦) 


PROCEDURE  MATCHPIC, 
LAPEL  767. 


MEAN). 


) 
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VAR 


STAT, P COUNTER  .  DNUM. PNUM,  DUL,  DLL  INTEGER, 

PROW. PCOL  0  255. 

DBLOCK. PJCBlOCK  PACKED  ARRAYCO  633  OF  DATAFILE, 

NFILE. pfile.  dfile  datafile* 

FNAME.  ORIGINAL,  P 1CNAME  DATAF, 


BEGIN 

767  WR1TELN('  ENTER  THE  CORRESPONDING  LETTER  FOR  THE  OPERATION  TO  BE  PERFORMED  - 
WRITELN. 

WR I TELN( '  M  A  ”  FOR  THE  *  AND  OPERATION  *  '), 

wRITELNt'  "  0  "  FOR  THE  *  OR  OPERATION  */), 

READLN* ANSWER ). 

IF  <  ANSWER  O  ANDCHAR)  and  (ANSWER  <>  DRCHAR)  THEN  GOTO  767. 

IF  ANSWER  *  ANDCHAR  THEN 
BEGIN 

WRITELN, URITELN, 

WR1TELN( 'YOU  CHOICE  OF  THE  AND  OPERATOR  REQUIRES  THAT  YOU 
WRITELN( 'ENTER  A  DECIMAL  NUMBER  FROM  0  TO  9  OF  THE'). 

WR  I TELN(  'STANDARD  DEVIATION  TO  BE  USED  EG  7  OR  3'). 

WR  I TELN (  '  THE  CHOICE  OF  0  WILL  INITIATE  AN  EXACT  AND  OPERATION  WITH  NO  THRESHOLD  WINDOW') 
READLN ( KONST ) . 

STAT  =*  ROUND  V  KONST  *  STDEVI  AT  I  ON  ); 

END 

ELSE 

BEGIN 

WRITELN. WRITELN. 

WRITELNt  'YOUR  SELECTION  -  THE  OR  OPERATOR  -  IS  NOW  EXECUTING'). 

END, 

F1LETITLE(FNAM£.  'Nl'>; 

FILETITLE(PICNAME,  PIC1'). 

FILETITLE(ORIGINAL.  Dl). 

PCOUNTER  »  1 i 

WHILE  PCOUNTER  <=  4  DO 

BEGIN 

IF  PCOUNTER  *  1  THEN 

WRITELN( 'PROCESSING  F 1 LE8LOCK  I'), 

RESET (PICNAME). 

RESET  ( OR  I G  I NAL  ) . 

FOR  PROW  =  0  TO  63  DO 
BEGIN 

GET {PI CNAME ) i 

P I C BLOCK! PROW 3  =  PICNAME^i 

GET (ORIGINAL); 

DBLOCK  [PROW  3  :»  ORIGINAL*4; 

PFILE  *  P I C BLOCK  C  PROW  3 . 

DFILE  =  DBLOCK L PROW  3 ; 

WITH  PFILE  DO 
BEGIN 

FOR  PCOL  -  13  TO  255  DO 
BEGIN 

IF  ANSWER  =•  ANDCHAR  THEN 
BEGIN 

WITH  DFILE  DO 
BEGIN 

DNUM  *  LI  NEC PCOL 3. 

DUL  *  DNUM  ♦  STAT; 

DLL  . *  DNUM  -  STAT; 

END, 

PNUM  =  LI  NEC PCOL 3. 

IF  ((DLL  O  PNUM)  AND  (PNUM  DUL ) >  THEN 
BEGIN 

IF  ANDTYPE  *  TRUE  THEN 
BEGIN 
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NF  ll_£  L  INt  l  PCOL  J 


END 

ELSE 

BEGIN 

NFILE  LINECPCOL3  =  PNUM, 

END 

END 

ELSE 

BEGIN 

NFILE  L1NEIPCOL3  -  1 5; 

END 

END  (*  ANSWER  *  ANDCHAR#) 

ELSE 

BEGIN 

WITH  DFILE  DO 
BEGIN 

DNUli  =  LINECPCOL3; 


END; 

PNUM  =  L I NE  t  PCOL  3  . 

IF  <  PNUM  =  15)  AND  (DNUM  =  15)  THEN 
BEGIN 

NFILE  LI  NEC PCOL 3  :  =  15. 

END 

ELSE 

BEGIN 

IF  DNUM  <  PNUM  THEN 
BEGIN 

NFILE  LINECPCDL3  .=  DNUM, 

END 

ELSE 

BEGIN 

NFILE  LINE  C  PCOL  3  .=  PNUM; 

END 

END. 

END. 

END. 

END. 

RESET (FNAME ) . 

SEEK (FNAME, PROW) . 

WR I  1 E ( FNAME,  NFILE); 

CL OSE< FNAME); 

END; 

CLOSE ( P I CN4ME  ) » 

CLOSE (ORIGINAL ) . 

IF  PCOUNTER  =  1  THEN 
BEGIN 

F  ILETI TLE ( FNAME.  ' N2 '  )  . 

FILETITLE(PICNAME.  'PIC*?')* 

F ILET ITl  E (ORIGINAL,  D2  '  )  , 

URITELN<  'FILEBLOCK  2  NOW  BEING  PROCESSED'), 

END 

ELSE 

IF  PCOUNTER  *  2  THEN 
BEGIN 

FILETITLE<  FNAME,  'N3'). 

P  ILETI TLE < P I CNAME,  PIC3). 

FILETITLE (ORIGINAL.  'D3'). 

WRI TELN( 'FILEBLOCK  3  NOW  BEING  PROCESSED'); 

END 

ELSE 

IF  PCOUNTER  «  3  THEN 
BEGIN 

FILETITLE<  FNAME  >  'N4'>. 

FILETITLE(PICNAME.  '  P  I C  4  '  ) , 

FlLETITLE( OR  I C I NAL.  ' D4 ' ) . 

WRI TELN( 'FILEBLOCK  4  NOW  BEING  PROCESSED’). 
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PCOUNTER 


PCQUNTER 


1. 


END. 

END. 


BEGIN  <*  MAINLINE#) 

WRITELN,  WRITELN.  WRITELN.  WRITELN,  WRITELN.  WRITELN.  WRITELN,  WRITELN. 

WRITELN. WRITELN. WRITELN. 

*R I  T£LN<  'WELCOME  TO  PROCESS  PICTURE  -  YOU  ARE  CURRENTLY  USING  THE  PASCAL  PROGRAM  . 
WRITELN. 

WRITELN!'  "ANDOR  PA"  (VERSION  MAY  19B2)'). 

WRITELN,  WRITELN.- 

WRITELN! 'ENSURE  ALL  PICTURE  FILES  ARE  PROPERLY  SETUP'). 

WRITELN! 'ON  DISK  BEFORE  PROCESSING  WITH  THIS  PROGRAM  IE.  FILES  MUST  EXIST')* 

WR 1TELN< 'FOR  PIC1  THRU  PIC4  AND  Dl  THRU  D4  (USE  FORTRAN  PROGRAMS  -  '). 

WRITELNt  "’VPIC"  AND  “VD"  TO  CREATE  THESE  FILES  IF  YOU  HAVE  NOT'). 

WRITELNt 'ALREADY  DONE  SO)  THE  PICTURE  FOR  WHICH  THE  PICTURE')* 

WRITELN<  MEAN  AND  STANDARD  DEVIATION  WILL  BE  COMPUTED  .  IF  DESIRED  BY  THE  USER.  MUST'). 
WRITELN(  'EXIST  IN  FILES  DAT  A 1  THRU  4  (VIA  "V4"  FORTRAN  PROGRAM)'); 

WRITELN, WRITELN. 

1111  WR1TELN(  'BESIDES  PERFORMING  "AND”  AND  '•OR”  OPERATIONS'), 

WRITELN! 'THIS  PROGRAM  CAN  ALSO  CALCULATE  THE  MEAN  AND  STANDARD  DEVIATION'). 

WRITELN! 'FOR  THE  PIXEL  VALUES  OF  THE  ENTIRE  PICTURE.  DO  YOU  WISH  TO')- 
WRITELN!  'HAVE  THESE  VALUES  CALCULATED?  Y/N  '). 

READLN! ANSWER); 

IF  ANSWER  »  YESCHAR  THEN 
BEGIN 

WRITELN.  WRITELN!'  ROGER,  USER  .  WAIT  ONE.  PLEASE'). 

WRITELN. WRITELN* 

THRESHOLD* 

WR  I TELN.  WR  I  TELN.  WR  I  TELN,  WR  1  TE LN; 

WRITELN!'  YUK  *  '  !  WHAT  A  PAIN  IN  THE  #£*7.  THAT  TURNED  OUT  TO  BE'"'), 

WR  I  TELN,  WR  I  TELN.  WR  I  TELN,  WR  I  TELN, 

FND 

ELSE 

IF  ANSWER  *  NOCHAR  THEN 
BEGIN 

WRITELN!  'ENTER  MEAN  AND  STANDARD  DEVIATION'). 

READLN(MEAN.  ST  DEV  I  AT  I  ON ) . 

WRITELN!  'MEAN  =  '.MEAN,  '  STANDARD  DEVIATION  *  ' .  ST DE V I  AT  I  ON ) . 

END 

ELSE 

BEGIN 

WRITELN!'  ERROR  !  RE-ENTER  VALUES' >* 

GOTO  1111. 

END. 

WRITELN,  UIR  ITELN.  WRIT  El  N, 

WRI7ELN<’D0  YOU  WANT  THE  OUTPUT  OF  THE  "AND"  OPERATION  TO  BE'). 

WR I TELN! 'B INARY  -  IE  GREYLEVEL  0  L  1 5  ONLY  -  (RECOMMEND  "N" )  Y/N')* 

READLN<  ANSWER  > . 

IF  ANSWER  *  YESCHAR  THEN 

ANDTYPE  .  -  TRUE 

ELSE 

ANDTYPE  «  FALSE. 

MATCHP IC. 

WR I TELN.  WRI TELN.  W» I TELN, 

WRITELN!'  PROCESSING  OF  ANDOR  IS  COMPLETE  '), 

WR  [  TELN. 

wRITFLNt  'PARAMETERS  USED  FOR  PROCESSING  WERE  AS  FOLLOWS  ')* 

-R  I  TELN.  WRITELN, 

-RITFcN!  ‘PICTURE  ME  AN  »  ',  MEAN)* 

-srfvN(  PICTURE  STANDARD  DEVIATION  *  STDEVI  AT  I  ON  ) , 

:>  -r.DT'FE  =  TRUE  ThFN 

-9I'U_N(  'OUTPUT  iriACE  VALUES  UL«E  SE  LtC  T  ED  TO  BE  BINARY,  JE  O  AND  15  ONLY'J. 

IF  ANDTYPE  =  FAt SE  THEN 

^RITE^N!  OUPUT  IMAGE  VALUES  ARE  Those  OF  PICTURE  STORED  IN  Dl  THRU  D4  '  ) ; 
wR I TELN!  'DEC IMAL  NUMBER  USED  FOR  THRESHOLD  WINDOW  »  '.  KONST), 

END 
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